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Abstract 



We have derived the most general chirally invariant Lagrangian Cq for the 
meson sector at order p^. The result provides an extension of the standard 
Gasser-Leutwyler Lagrangian £4 to one higher order, including as well all 
the odd intrinsic parity terms in the Lagrangian. The most difficult part of 
the derivation was developing a systematic strategy so as to get all of the 
independent terms and eliminate the redundant ones in an efficient way. The 
'equation of motion' terms, which are redundant in the sense that they can be 
transformed away via field transformations, are separated out explicitly. The 
resulting Lagrangian has been separated into groupings of terms contributing 
to increasingly more complicated processes, so that one does not have to deal 
with the full result when calculating contributions to simple processes. 
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I. INTRODUCTION 



Ever since the early days of current algebra and the PC AC hypothesis (for an overview 
see, e.g., Refs. approximate chiral symmetry and its application in terms of effective 

Lagrangians have been cornerstones of the description of the low-energy interactions 
of hadrons. In our present understanding, chiral symmetry is interpreted in terms of QCD, 
the SU{3) gauge theory of the strong interaction involving quarks and gluons (see, e.g., 
Refs. In the limit of massless u, d and s quarks, the QCD Hamiltonian exhibits a 

global SU (3) L X SU (3)/j symmetry which is assumed to be spontaneously broken to SU (3)y, 
giving rise to 8 massless Goldstone bosons. The experimentally observed small masses of 
the pseudoscalar octet (vr, rj) then originate from an explicit symmetry breaking due to 
the finite quark masses (see Refs. [0,|n|] and references therein). 

Extending a method originally proposed by Weinberg for the analysis of S-matrix ele- 
ments [^], Gasser and Leutwyler |T3,n] developed a technique which allows an expansion 
of QCD Green's functions of quark currents in terms of momenta and quark masses. Their 
procedure, now known as chiral perturbation theory (ChPT), makes use of an effective La- 
grangian approach for the interaction between the Goldstone bosons vr, K^rj, and has been 
applied to a wide variety of processes, including interactions with external electromagnetic 
and electroweak probes (for pedagogical introductions and recent reviews, see Refs. [|T3]-^). 
The lowest-order Lagrangian corresponds to a nonlinear a model containing two free pa- 
rameters, namely, the pion decay constant and the scalar quark condensate. At 0{p^) the 
most general Lagrangian which is consistent with chiral symmetry, parity, and charge conju- 
gation invariance contains 10 structures. The renormalized coefficients Li . . . Liq have been 
empirically determined by fitting experimental data |13,l3 (fo^' ^ recent analysis see Ref. 
and references therein). Theoretical predictions for the chiral coefficients have also been 



obtained by several authors using various techniques (see, e.g., Refs. |r3| , |23| -p8|]). 

Without external fields (i.e., pure QCD) or including electromagnetic processes only, the 
effective Lagrangian of Gasser and Leutwyler has an additional symmetry resulting in the 
property that it contains interaction terms involving exclusively an even number of Goldstone 
bosons pP|J5(][] . Such interaction terms are sometimes also referred to as being of normal or 
even intrinsic parity. In Ref. [|29| Witten discussed how to remove this symmetry which is 
not a symmetry of nature (for example, 7r° -j'j,K~^K'' tt'^it^tt^, etc.). He essentially 



rederived the Wess-Zumino anomalous effective action describing the chiral anomaly |^ 



The corresponding Lagrangian, which is of O(p^), cannot be written as a standard local 
effective Lagrangian in terms of the chiral matrix U but can be expressed directly in terms 
of the boson fields. In particular, by construction it contains interaction terms with an odd 
number of Goldstone bosons (odd intrinsic parity). It was subsequently shown by several 
authors that quantum corrections to the Wess-Zumino classical action do not renormalize 
the coefficient of the O(p^) Wess-Zumino term p2|-|5^. Furthermore, the one-loop counter 
terms lead to conventional chirally invariant structures at 0{p^) |33|-p^. For a review of 



chiral perturbation theory involving the odd intrinsic parity sector we refer the reader to 
Ref. g. 

Chiral perturbation theory to O(p^) has become a well developed effective representation 
of QCD at low energies in terms of Goldstone bosons, with many applications. However, 
there exist cases where one-loop calculations to 0{p'^) do not appear to lead to satisfactory 
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agreement with experimental data (see for example Refs. ||37| , |38[] ). As consistent two-loop 
calculations in the even intrinsic parity sector for particular processes have started to become 
available it is clearly important and timely to extend the most general chiral effective 
Lagrangian to 0{p^). Moreover, it is now also possible to calculate the O(p^) coefficients 
from phenomenological models pO|] . 

The purpose of this work is to derive the most general structures of the chiral Lagrangian 
at 0{p^) including the even, as well as the odd, intrinsic parity sector. The result provides an 
extension in analogous form of the standard Gasser-Leutwyler O(p^) Lagrangian to 0{p^). 

In the even intrinsic parity sector we find 111 independent terms. This part of the 
0{p^) Lagrangian will be of relevance for processes which vanish at O(p^), so that O(p^) 
is the leading order and the 0{p^) terms provide the leading-order correction. A two-loop 
calculation of such a process, namely 77 tt^tt^, has recently been reported in Ref. (see 
also 10). It is a feature of this process, and similarly of the decay i] — >• 7r°77 |3^, that the 
OIp"^) contribution is exclusively generated by one-loop diagrams. According to the power 
counting scheme [0, at 0{p^) contributions result from two-loop diagrams with vertices 



from the 0{p'^) Lagrangian, one-loop diagrams with one vertex being 0{p'^) and the other 
0{p'^) and, finally, from tree-level diagrams from the 0{p^) effective Lagrangian. 

In the odd intrinsic parity sector we find 32 independent structures of 0{p^). Our result 



differs from previous determinations of these structures |34J3q] which in turn are not in 



agreement with each other. We will provide a detailed comment on these differences. This 



sector has already been discussed in the hterature in quite some detail as it provides 
leading-order corrections to processes originating in the Wess-Zumino term. 

At present it appears extremely unlikely that all renormalized coefficients at 0{p^) may 
be determined empirically. However, there is the hope that for simple processes the much 
smaller subgroups of relevant coefficients can be determined by one experiment and then 
provide predictive power for a related process. Furthermore, theoretical techniques which 
resulted in predictions of the 0{p'^) coefficients might be extended to predict the coefficients 
at 0{p^). For the corresponding discussion of techniques applying to the odd intrinsic parity 
sector see Ref. |3^. A derivation of the 0{p^) even intrinsic parity sector using the Nambu- 



Jona-Lasinio model will be provided in Ref. ||40|| . 

Our main emphasis in this paper will be on developing a systematic procedure for the 
derivation of the 0{p^) Lagrangian so that the reader can be confident that the large number 
of terms we have found are indeed independent as well sufficient to describe the most general 
effective chiral Lagrangian at 0{p^). The derivation involves various technical points, such 
as total-derivative arguments, trace-relations, use of the classical equation of motion in 
terms of field transformations, and special relations involving the completely antisymmetric 
tensor in four dimensions. We group the final structures according to the minimal number 
of Goldstone boson fields, assuming for the purpose of organizing the terms, a coupling to an 
external electromagnetic field. Our list will allow the reader to quickly identify the relevant 
terms at 0{p^) needed for particular physical processes. The structures will, of course, be 
given in their full generality without any assumption concerning the nature of the external 
fields. 

Our work is organized as follows. In Sec. y we discuss some of the general foundations 
of chiral perturbation theory, define the basic building blocks for the Lagrangian, discuss 
the strategy for obtaining a complete and independent set of terms, and describe a number 
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of simplifications, symmetries and relations which will be used to reduce the result. In Sec. 
PI we describe the details of the explicit construction of the Lagrangian. Sec. IV contains 
a discussion of the equation of motion terms, of some simplifications, and of a scheme of 
organization which groups the most useful terms together and a comparison with some 
previous works. The final results are given in Tables ^ - |VI1| . Sec. |V| gives a brief summary. 
Details of the derivation of trace relations used in the simplifications, a table of the equation 
of motion terms, and some tables showing relations among terms which are not independent 
are relegated to the appendices. 

The reader primarily interested in using the results should read Sec. || to understand 
the notation and the general strategy and Sees. |rV| and |V| for a listing and discussion of the 
results. 



II. CONSTRUCTION OF THE MOST GENERAL LAGRANGIAN Cq - 
INTRODUCTION AND BASIC APPROACH 

A. Introduction 

Our aim in this section is to lay the groundwork for the construction of C^, the most 
general chirally invariant Lagrangian at order p^. This Lagrangian will provide an extension 
to the next higher order of the standard Gasser-Leutwyler £4, as well as the odd intrinsic 
parity Lagrangian. Like £4 it will be expressed in terms of a set of basic building blocks, 
arranged in traces or products of traces such that the result is chirally invariant. 

The main problem with such a construction is that it is far too easy to think of terms 
satisfying the necessary criteria. One quickly obtains many more terms than necessary, 
many that are not independent, though often not obviously so, and many expressed in a 
more complicated fashion than necessary. 

Thus the most difficult task in such a construction is to develop a systematic strategy 
which allows one to get in an efficient manner all of the independent terms while at the 
same time eliminating, or preferably never generating, those terms which are redundant. 
Most of this section will be devoted to developing the various ingredients of this strategy. 
Then in the following major section the strategy will be applied to derive the most general 
Cq. Some of this will be pedagogical, as it is important for the reader to be able to follow 
and understand the steps of the derivation and thus be convinced in the end that the result 
obtained is, in fact, complete and correct, something which is not necessarily possible when 
only a table of results is given. 

To our knowledge there have been no previous attempts to generate the complete and 



most general C^, though there have been several attempts p^j36[1 to derive the set of odd 
intrinsic parity terms, i.e., those involving the completely antisymmetric tensor eap-ys- In 
all of these attempts it appears that terms were missed, redundant terms were included, or 
both. 
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B. Chiral symmetry of the QCD Lagrangian 



Before discussing the construction of the effective chiral Lagrangian we will shortly review 
the chiral symmetry of the QCD Lagrangian coupled to external fields. This discussion is 
relevant for the identification of the transformation properties of the external sources under 
the group, parity, and charge conjugation. 

Let us consider the QCD LagrangianQ 



1 



f 



(1) 



In Eq. (|1]) qf denotes a quark field of flavor / with current quark mass rrif. The covariant 
derivative is defined as 



(2) 



where, for simplicity, we do not exhibit the fact that the quark fields have three color indices. 
In Eq. (Q) g denotes the strong coupling constant. A" are the usual Gell-Mann matrices, and 
represents a gluon field with color index a which may take the values 1 to 8. Summation 
over repeated color indices is assumed. The gluon field strength tensor is defined as 



(3) 



where fahc are the SU (3) structure constants. We do not discuss the 6'-term which induces 
P, T and CP violations in the strong interactions 0. 

In the following we will restrict ourselves to three flavors, m, d, s, and consider the limit 
mu,md,ms ^ 0. In this limit Eq. (|I|) reduces to 



^'qcd= E {^qLJrD,qLJ + ^qRJrD,qHJ)--C%Gl^'', (4) 

f=u,d,s 

where the left-handed and right-handed quark fields are defined as 

QL = ^{i-75)q, qR = ^{^ + i5)q- 

The Lagrangian of Eq. is invariant under the following global transformations of qi and 
qR, respectively. 



d'L 



Ur 






(5) 



We omit the gauge fixing term and the Faddeev-Popov ghost term, since they are not relevant 
for the discussion of chiral symmetry. 
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where Ul and Ur are independent ^7(3) matrices. This invariance, in principle, gives rise 
to 18 conserved currents. However, due to quantum effects, the axial U{1)a current is 
not conserved. In the following we will only be concerned with the analysis of the G = 
SU{3)l X SU{3)ji symmetry of Eq. (Q). It is generally accepted that this symmetry is 
spontaneously broken to SU{3)v giving rise to 8 Goldstone bosons. The finite quark masses 
in Eq. (|^) then give rise to finite - but in comparison with other hadrons small - masses of 
the Goldstone bosons. 

In order to systematically study the consequences of chiral symmetry and its breaking 
through the quark masses, we follow the technique of Gasser and Leutwyler |r3,n] and 
introduce external c-number fields v^{x), a^{x), s{x) and p{x) into the Lagrangian, 

^ = ^QCD + ^ext = ^QCD + Ql^i'^t^ + 75^)^ " q{s - ^TsP)^- (6) 

The external fields are color neutral, hermitian 3x3 matrices, where the matrix character, 
with respect to the (suppressed) flavor indices u, d, s of the quark fields, is 

A'* A'* A'* A** 

Of course, the three flavor QCD Lagrangian is recovered by setting = = p = and 
s = diag^niu, md, mg). 

Requiring the total Lagrangian of Eq. (j^) to be invariant under P, C and T leads to con- 
straints on the transformation behavior of the external fields. Under a parity transformation 
the quark fields transform a^ 

qf{x,t)^^\f{-x,t). (8) 

The requirement of invariance under a parity transformation, 

L{x,t) ^ C{-x,t), (9) 

results in the following transformation properties of the external fields 

v^, -a^, s ^ s, p ^ -p, (10) 

where it is understood that the arguments change from {x,t) to {—x,t). Under charge 
conjugation the quark fields transform as 

c c \ 

QaJ ^ Caf3q(3J, QaJ ~^ ^QfijCp^. (11) 

In Eq. ( pUD the subscripts a and are Dirac spinor indices, C = 27^7° is the usual charge 
conjugation matrix in the convention of Ref. and / refers to flavor. Using Eq. ( pA]) 



it is straightforward to show that invariance of Cext under charge conjugation requires the 
transformation properties 



^We suppress color indices since we are considering color neutral external sources. 
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C T C T C T T f-\n\ 

v^^-v^, a^^a^, s,p^s,p, (12) 

where the transposition refers to the flavor space. 

Finally, we discuss the properties of C under the group G. To that end we rewrite Eq. 
(I) in terms of ql and qr, 

^ = ^QCD + QlI^L^ql + QR'^^Rt.qR + Qhis - ip)qR + ^/?(s + ip)qL, (13) 

with Rfj^ = + and = — a^. We then promote the global symmetry to a local 
symmetry. Eq. (p!3|) remains invariant under qr — > VR{x)qR and q^ ^^(x)^^, provided the 
external fields transform as 

^ VrR.V^ + iVrO^V^, ^ VlL.vI + iVlO^vI 
s — ip ^ Vl{s — ip)VR, s + ip Vr{s + ip)Vl. (14) 



C. Basic Building Blocks 

The Lagrangian is constructed from the same basic ingredients as £4, namely the 
Goldstone boson fields, external gauge fields, scalar and pseudoscalar external sources, and 
their derivatives. The 8 Goldstone bosons arising from the spontaneous symmetry breaking 
are collected in a SU (3) matrix 



U{x) 



exp 




(15) 



with 



-71 



75^ 



(16) 



V2K0 J 



and Fq the pseudoscalar meson decay constant in the chiral limit WM. The matrix U 
transforms linearly under the group G = SU{3)l X SU{3)r : U ^ U' = VrUVI 



Furthermore, under charge conjugation and parity the Goldstone bosons transform as G : 
(;/)—> 0^ and P : (p{x,t) ~(p{~x,t), or equivalently G : U ^ U'^ and P : U{x,t) — >■ 
[/t(-f, t). 

We define the field strength tensors associated with the external gauge fields R^ and 

as 



= (^ti^u — duR^i — i[Rfi, Ru], 



(17) 



Note that and F^^, are hermitian, traceless, and gauge invariant, i.e., under the group 



they transform as ^RFj^^VR and ^-^/fj/^- Furthermore, we follow Ref. [14 
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and introduce the linear combination x = 2-Bo('S + zp), with x ^ yRX^li where Bq is related 
to the vacuum expectation value < 0|gg|0 >. 

The effective Lagrangian is constructed in terms of [/, ?7^, x, x^ ^ind the field strength 
tensors F^, as well as covariant derivatives of these objects. These covariant derivatives 
involve the gauge fields and and transform in the same way under the group as the 
quantities they act upon. Given the transformation properties of and of Eq. ( p!4D we 
define the covariant derivatives as 



A ^ VrAvI 
B ^ VlBVI 
C ^ VrCV}, 
D ^ VlDVI 
E 



D^A = d^A-iR^A + iAL^, 
D^B = d^B + iBR,, - iL^,B, 
D^C = d^C-iR^C + iCR^, 



D,E 



d,E. 



(18) 



Note that we use the same symbol for the covariant derivative, independent of the 
transformation property of the object it acts upon. The advantage of this convention is that 
a chain rule analogous to ordinary derivatives holds. Given the product Z = XY where 
X, Y, Z have, according to Eq. (0), well-defined but not necessarily the same transformation 
behavior, the chain rule applies 



Df,Z 



D^{XY) = {D^X)Y + X{D^Y), 



(19) 



which is straightforward to verify using the definitions of Eq. 

This chain rule is valuable as an intermediate step in a number of the derivations of 
various relations. In essentially all cases however the final results will be expressed solely in 
terms of quantities transforming as U and their covariant derivatives and thus require only 
the covariant derivative defined in the first line of the definitions of Eq. ([181) . 

It will be very useful in the subsequent derivations to have all of the building blocks 
transform in the same way and also to be able to handle the external field terms Xi -^/l? -^/ft/ 
in the same way. To that end we define 



R 



uF^r 



(20) 



We then let x'^'^ = X ^i^nd let x'^'^, with jj, ^ stand for G^'^ or H^'^ . Thus we can consider 
X^^ with /i, V unspecified to stand for, collectively, the set (x, G^" , H 



With these definitions, we have only two basic building blocks U, x^^ and their adjoints 
and covariant derivatives acting on them. Note that by virtue of the chain rule, we do not 
need derivatives acting on products of these basic terms. All building blocks then transform 
as U (or [/"I"). In terms of the momentum expansion U is of order 1, x'^'^ is of order and 



each covariant derivative is of order p |T^ . 

Finally, in order to construct terms which are invariant, we must put these basic pieces 
together so as to get quantities which transform as . . . V^. Thus define for any A trans- 
forming as VrAVI, i.e. as U, the hermitian (or antihermitian) combinations 
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[A]^^-iAU^±UA^). 



(21) 



We can then take as the basic building blocks [A]± with A taken as x'^^ or D^U, or as 
some number of covariant derivatives acting on x'^'^ or D^U. The trace of a string of such 
quantities, or the product of such traces, will then be invariant under the group. Observe 
also that by defining the building blocks this way we include each operator and its adjoint. 
By using building blocks transforming as Vr . . . we also solve the problem of determining 
whether or not extra f/'s and f/^'s are necessary. In the usual approach, starting with a 
string of quantities transforming as Vn . . .Vl or Vl ■ ■ ■ V^, one can always insert extra t/'s 
and W^s in the string, replace some terms by their adjoints, and get another invariant string. 
One then must check in each case whether the resulting string is in fact independent of the 
original one. However by using the forms of Eq. ( |2TD the only possible combination which 
can be inserted is UW which is unity and of course unnecessary. Observe also that if we 
were to insert UW between each term and after the last term in a trace of a string of terms 
of the form of Eq. (|lD and use the relation W{AW ± UA^)U = WA ± A^U we obtain a 
trace of the same form as the original, but constructed of building blocks transforming as 
Vl ■ ■ - Vl which we could have defined in a fashion analogous to Eq. ([2l|) . This shows that 
such terms are not independent and that terms of the form of Eq. (^) are sufficient. 

To summarize, we have seen that the basic building blocks necessary to construct the 
most general chirally invariant Lagrangian are of the form {Al\±, where A is -D^t/, x^'^ 
(multiple) covariant derivatives of such quantities. The problem now, to be addressed in 
the next section, is how to put these building blocks together in an efficient way to get all 
possible terms in the Lagrangian. 

D. Strategy 

It is possible to generate a very large number of terms using just the building blocks 
defined in the previous section. The most difficult part of obtaining is developing a 
strategy which obtains all of the independent terms without generating a lot of extraneous 
terms which have to be eliminated by hand. 

To do this it is convenient to first define a hierarchy of terms. At each level we will then 
find the most general set of terms. First however, various relations can be used to eliminate 
some terms in favor of those lower down in the hierarchy. Since at each level we always find 
the most general set, it is not necessary to actually work out what these (often extremely 
complicated) relations are. To eliminate a term one must only show that there exists a 
relation expressing it in terms of quantities being kept at the same level and others lower 
down in the hierarchy. 

We thus order the various classes of terms as follows: 

I. Those terms with six -D^'s. 
II. Those terms with four -D^'s and one x^^ ■ 
III. Those terms with two -D^'s and two x^'^'s- 



9 



IV. Those terms with three x^'"s. 

At each major level terms which involve a single trace are considered higher up than those 
with multiple traces. Likewise the ordering ensures that the addition of an F^'^ always 
results in a term lower in the hierarchy. Since the result must be a Lorentz scalar and since 
the only available tensors ^ind eap-ys have an even number of indices there are no terms 
with an odd number of D^'s. 

For each level one first determines the form of the possible structures. For this it often is 
easier to start initially with the original building blocks U, x^'^ ^"^^ their covariant derivatives. 
A number of simplifications and tricks to be discussed in the next section are then applied 
to reduce the number of possible terms and ensure that each term satisfies parity and charge 
conjugation. The remaining structures are then written in terms of the final building blocks 
{A\±. The result is a set of general forms. Explicit results can then be read off directly, 
usually by simply considering all possible permutations of the indices, and perhaps of the 
order of the terms, and all possible ways single traces can be broken up into multiple traces. 

Terms involving D^D^U are of particular interest as this factor appears in the classical 
equation of motion and terms proportional to the equation of motion can, as we shall see 
in a later section (see also Ref. [Q), be transformed away by a transformation on, or a 
redefinition of, the fields. Thus our strategy will be to extract in so far as it is possible such 
terms, which will be referred to as 'equation of motion terms'. At the end we can make such 
terms proportional to the full C2 equation of motion by simply adding in the rest of the 
terms in the equation of motion, as these terms will come from a lower level in the hierarchy. 



E. Simplifications-General Results 

1. Total derivative arguments 

It is convenient to make use of the fact that a total derivative in the Lagrangian density 
does not change the equation of motion. This 'total derivative argument' can be applied 
as follows. For any pair of operators A and B which under the group either transform 
as A = VrAVI and B' = VlBV^ or as A' = VrAV^ and B' = VrBV^ the product AB 
transforms as VrABVr and the chain rule and definitions of the covariant derivatives yield 

D,{AB) = {D^A)B + A{D^B) = d^{AB) - i[R^, AB], (22) 

and thus using the fact that the trace of a commutator vanishes, 

Tr {{D^A)B) + Tr {A{D^B)) = d^Tr (AB) . (23) 

The same result is obtained for combinations AB transforming as VlABVI- 

Clearly this result generalizes to products of terms of the form Bi . . . Bn, and one gets 

Tr {{D^B,)B, . . . 5„) + Tr {B,{D^,B,)B; . . . 

+ . . . + Tr (Si . . . (D^B^)) = d, (Tr (B, . . . B^)) . (24) 

It is also straightforward to extend the 'total derivative argument' to products of traces. 
As an illustration, we find applying Eqs. (El) and (|T^), 
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Tr {iD^A,)A2) Tr + Tr (A.iD^A^)) Tr 

+Tr Tr {{D^A^)A^) + Tr {A^A^) Tr {A^{D,A^)) 

= d,{Tr{A,A2)Tr{A^A^)). (25) 

As a consequence of these arguments we may move covariant derivatives from one fac- 
tor to another either within a single trace or across muhiple traces in accord with Eqs. 
(p3|) , (p^) , (p5D ■ Terms involving total derivatives will result, but these will not contribute to 
the equation of motion and thus can be dropped from the Lagrangian. 



2. Symmetrization of multiple covariant derivatives 
For any operator A transforming as U we have, using the definition of covariant deriva- 



tives Eq. (p!8| ) and of the field strength tensors Eq. (|T1 

{D,D, - D,D,)A = zAF;:, - iF^^A. (26) 

Any term of the form D^Di^A can be written as a term symmetric and one antisymmetric 
m ^ ^ V. Thus Eq. (^) implies that at any given level in the hierarchy D^D^A can always 
be assumed to be symmetric in its indices, since the antisymmetric part can be expressed in 
terms of the field strength tensors and thus contributes only at a lower level in the hierarchy. 
Clearly this can be generalized to more than two covariant derivatives. Hence we will always 
assume that a multiple covariant derivative, D^^ . . . D^^A, is symmetric in the interchange 
of any of its indices. Formally this means that we will always interpret D^D^jA as 

D^D^A ^ ]^{D^,D, + D,D^,)A, (27) 
with an analogous definition when there are more than two covariant derivatives. 



3. Index exchange 

In keeping with the strategy of extracting as many 'equation of motion' factors as possi- 
ble, one would like to arrange things so that covariant derivatives with summed indices, e.g. 
DaD", act on the same factor U. For example consider Tr {{DaDf^A){D°'DjB)C), where 
A, B, C do not contain any covariant derivatives. Using the results of Sees. |11E 1| and [11 E 2 
one can move the off of A and on to B and C and then move the Da off of B and on to 
A and C . The result is 

Tr {{DaDpA){D^D^B)C) = Tr {{DaD^A){DpD^B)C) 

+terms with multiple D's on only one of A,B,C 
+total derivatives, (28) 

where we have dropped the terms arising from the commutation of the D's at various stages 
which can be expressed via terms lower in the hierarchy using Eq. (pG]). 
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4- Building blocks involving multiple covariant derivatives 

Consider a term of the form [DaDpU]^. By starting with the trivial identity = 
DaDpiUW) and using Eq. (|19D one finds 

2[D^DfsU]+ = -{D^U){Dfs0 - {DpU){D^U)^ (29) 

which means that [DaDpU]^ can always be expressed in terms of quantities involving at 
most a single covariant derivative acting on U . This can be generalized, starting with 
= Dfj_^ . . . D^^{UW), to show that [D^j . . . can always be expressed as a sum of 

products of [A]± where A has at most one less covariant derivative. Symbolically, [D^U]^ ~ 
{D'U){D''-'Uy ~ [D'U]±[D''-'U]± for all {I < i < n-1). By iteration, one then shows that 
all terms of the form [D^-^ . . .D^^U]^ can be expressed via the combinations with relative 
negative sign [A]_ and the right (or left) hand side of Eq. (p9D. Note also that the simplest 
case [Dfj_U]+ = 0, so that the right hand side of Eq. can also be expressed solely in 
terms of [A]^ as well. This means in practice that we never need to include factors of the 
form [D,,...D^„UU. 

5. Symmetries, parity and charge conjugation 

Each term in the final Lagrangian must be hermitian and must be invariant under parity 
and charge conjugation. Hermiticity is normally trivial as all terms come out to be either 
hermitian or antihermitian once parity and charge conjugation invariance is imposed. The 
strategy for ensuring P and C invariance is to take each candidate term and add to it the 
parity transform of the term, thus getting a result which is trivially P invariant. Similarly 
one then adds to this result its C transform, and so generates a C and P invariant term. 

Consider first parity. We will be interested in candidate terms for the Lagrangian which 
are of the form of a trace of a string of factors of the canonical form [A]±, or alternatively 
a product of such traces, where A can be Dfj^U or x'^'' or their (perhaps multiple) covariant 
derivatives. Table | shows the transformation properties of the various factors. We see that 
under parity all allowed A's transform as A —>■ (—1)^74''" , where in addition all Lorentz indices 
in A^ have been raised (or lowered) with respect to those in A, and x —x. The quantity 
p is an 'intrinsic parity' which is 1 for H^^ and its covariant derivatives and otherwise. In 
most cases Lorentz indices from one factor will be contracted with those of another factor, 
so that the raising or lowering of indices under P compensates in the two factors, and can 
be ignored. The exception to this is the case when the indices are contracted not with 
another term of the set of A's but with the antisymmetric tensor eap-ys- Parity does not 
affect this tensor, so to get a Lorentz invariant quantity one must raise or lower its indices 
by hand. This introduces a minus sign since e^^^^ = — e"^i"'. This can be summarized neatly 
by including a (—1)"^ in the parity transform of a complete term, where e counts the number 
of e^/j^^'s in the term. 

Using these preliminaries, we see that under parity [Ai\± -^{—1YU'^[A\±U . Under a 
trace the U'^ and U collapse to 1, using the cyclic properties of a trace. Thus we arrive at 
the following operational method of making a candidate term in the Lagrangian invariant 
under parity. Simply multiply the term by the factor 1 + (— 1)'^+^+^ where P is the sum of 
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the intrinsic parities p of the various terms, e is the number of Ca/j-y^'s in the term, and where 
the (—1)* is determined from the product of signs coming from the individual factors in the 
term, counting + for [A]_^. and — for Clearly this holds for a single trace or a product 

of several traces. 

Note also that by extracting this factor, we make it easy to see which terms are not 
allowed by parity conservation, namely those for which (— 1)*+-^+^ is —1. 

Consider now charge conjugation and proceed in exactly the same way. From Table 
I it is easy to see that the necessary A's transform under C as {—lyA^ where c is an 
'intrinsic charge conjugation quantum number' which is 1 for G'^'^ and its covariant deriva- 
tives and otherwise and T is the transpose. We then have under charge conjugation 
[A]± — > {—1Y{W[A]±U)'^. Using the fact that the trace of a quantity is equal to the trace 
of its transpose, we then find that under C the trace of a string of basic factors [A]± goes 
simply to (—1)'" times the trace with the basic factors in reverse order, where here C is 
simply the sum of the intrinsic charge conjugation quantum numbers of the various terms. 
Clearly in the case of multiple traces the factors are reversed in each trace individually. 

Thus to summarize, the operational way to get a term in the Lagrangian to be invariant 
under charge conjugation is to add to it a term (— 1)*" times the term with the major factors 
reversed in each trace. If all of the traces have at most two factors, the cyclic property of a 
trace implies that the second term has the same structure as the first, and that effectively 
one is simply multiplying by the overall factor (1 + (— 1)*^). 

As noted ensuring hermiticity is relatively simple since the basic building blocks satisfy 
[A]l- = ±[A]±. Thus the adjoint of a string of such terms is just (—1)^ times the string in 
reversed order, where s is calculated as in the parity discussion above. But, by virtue of C 
invariance, reversing the order produces an overall phase (—1)*^, so that taking the adjoint 
of a term simply generates an overall factor (—1)^+'^ = (—1)^+*^+'^. Thus to get a hermitian 
term we need only multiply by a factor of i when this factor is —1. This happens for terms 
with an even (odd) number of G^'^^s and iJ^'^'s (or their covariant derivatives) combined and 
one (no) ea/^^ys factor. We will do this only in the final listing of results. 



6. Trace relations 

In the construction of £4 a relation between SU (3) generators was used to eliminate one 
term . That relation can be generalized and used to eliminate a number of possible terms 



of order p^. A proof of the relation and some discussion is given in Appendix |^. For the 
present purposes we need only the result in the following form. For any set of 3 x 3 matrices 
Bi, i = 1, . . . ,n with n > 4, we can express 

Tr{B,...B^) (30) 

All perm. 

as a sum of products of traces, each trace containing no more than n — 1 of the Bi. Since 
we are interested only in terms invariant under the group, and since we take the Bi in all 
permutations, we need Bi to transform as VnBiV^. Thus the Bi of interest will be the basic 
building blocks [A]±. Since a product of traces is always lower down in the hierarchy than 
a single trace, this relation allows us to discard one single trace term from each set of terms 
made up of all permutations of a particular group of building blocks. 
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7. Epsilon relations 



There is also a set of relations which relate various terms involving the tensor eajj^ys, 
and allow one to eliminate some in favor of others [^]. These relations originate from the 



observation that a tensor antisymmetric in five Lorentz indices must be zero since there are 
only four possible different indices. Such a tensor is 

Now consider a further tensor Qap^yprr) which will actually be the trace, or product of 
traces, of the basic building blocks [A\± and contract this tensor with that of Eq. (^) in all 
possible ways. This leads to the following six equations: 

( Qa 'ypTTj ~^ Qa-y p^-^ Qa-yp ~l" Qaypr Qayprrj ^ 0, (32) 

(+Qa 'yprri Qya + Q^ap Q-yapr ^ + Q-yaprrj ^^^'^ = 0, (33) 

{,~\~Qayp Qjap + Qjpa rrj Qypra ^ + Qyprari ^^^^ = 0, (35) 

(—D -\- (1 ^ — O ° — (1 °')p^P^v — 

(_i_r) °' — (1 "-i-n ° — D "-l-D ° \f:-)p-rv — n /'Q7^ 

Observe the symmetry of the Q's across the diagonal and the fact that Eq. (^) is just the 
sum of the first five. Thus in general there are five independent equations which generate for 
each specific Qap-yprr] five relations among the terms involving e"^'^^ . However if Qapyprri has 
some symmetry in its indices, not all of these five relations are independent. For example 
if Qapyprr] is Symmetric in the first two indices then Eqs. (0), (|35|) , and (^) are identities 
and Eqs. ( |5^ ) and (^3p are not independent and so there is only one independent 'epsilon 
relation'. 



III. CONSTRUCTION OF THE MOST GENERAL LAGRANGIAN Cq - 

EXPLICIT CALCULATION 

Having completed these preliminaries we now proceed to the explicit calculation of C^. 
As noted earlier, the strategy will be to start with the highest level in our hierarchy and 
work down. At each level we first itemize the types of terms possible. At this stage it is 
useful to adopt a schematic notation. Thus an expression like {D'^U){DU) will mean the 
class of terms which have a {D^D^U) or its adjoint and a DaU or its adjoint, in any order, 
and with any contraction of indices which lead to a chirally invariant term when the trace is 
taken. We then use the results of Sec. [11 E| to eliminate as many classes of terms as possible 
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and to force invariance under parity and charge conjugation. The results are a set of general 
expressions which satisfy the various constraints. One can then read directly from these 
expressions all of the allowed terms by simply taking all possible index contractions, and in 
some cases all possible orders. 

At a given level in the hierarchy the manipulations required to reduce terms to lower 
levels and to ensure P and C invariance are essentially the same for multiple trace terms 
as for single trace ones. Hence we will generally work, somewhat symbolically, with just a 
string of factors and only at the end take the single trace and all possible multiple traces. 

We will always attempt to extract a factor [D^D'^U]- as the terms containing this fac- 
tor can be made proportional to the classical equation of motion, and as discussed below, 
transformed away via a transformation of the fields. 



A. Terms with six -D^'s and no x^'^'s 

Terms in this class will have six covariant derivatives acting singly or in multiples on 
U. Since there are six Lorentz indices, two must be contracted and the other four can be 
contracted pairwise in two pairs, or contracted with the indices of the antisymmetric tensor 
^a/B'yS- Using the total derivative arguments of Sec. [11 E 1| the covariant derivatives can always 
be moved back and forth, so as to get a factor (D^D^U) which will be pulled to the first of 
the string of terms which eventually will become an argument of a trace. Either = u in 
this factor or, if not, no other multiple derivative has a contracted index in it, since if it did 
we could move all but the two contracted D's onto another term and pull the contracted 
pair out as the first term instead. Note that by virtue of Sec. |11 E 2| all multiple derivatives 



can be taken as symmetric in their indices. This leads to structures of the following types: 

1) {D^DMKD^U), 

2) {D^D,U){D'U){DU)U, 

3) {D^DM){D^U){D^U)U, 

4) {D^DM){D'U){DU){DU), 

5) {D^D,U){DU){DU){DU){DU)U. (38) 

To get an expression transforming as V/j . . . so that the trace will be invariant we 
use Sec. [II E4| to write each of the structures in terms of products of factors [A\±. Each 



combination + , again by the results of Sec. [11E4[ , can be expressed via terms with one 



fewer covariant derivative, and hence via structures lower down in the list of Eq. (^), and 
so does not need to be kept. However at the bottom of the list [D^Di,U]+ is expressed in 
terms of {DU){DU) and so generates a sixth structure involving six [D[/]_'s. 

Requiring parity invariance introduces an overall factor 1 + (—1)^+^+^ in accord with 
Sec. [1IE5[ . The intrinsic parities p are all zero so we find that the parity conserving terms 
must have an even number of [A]^ factors and no eajj-ys or an odd number with an 
Finally charge conjugation invariance simply adds in a term with the [A]_ factors in reverse 
order. 

Thus we obtain terms of the following forms: 
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1) [D'UUD'U]^, 

2) [D'UUD'UUDUU^p,s, 

3) [D^UUD^UUD^UU^p^s, 

4) [D^U]4D^U]4DU]4DU]^, 

5) [D^UUDUUDUUDUUDUUap-yS, 

6) [DU]4DU]4DU]4DU]4DU]4DU]^, (39) 

where in each case we must add a piece with all major factors in reverse order to ensure 
charge conjugation invariance. We must also take all possible starting orders, which in 
this case, once C is taken into account, is relevant only for 4) and generates two different 
structures there. 

Note that all manipulations so far hold for single or multiple traces, except that for 
multiple traces we reverse the terms in each trace separately. 

Finally we put in indices and the order reversal explicitly to obtain the most general result 
for single trace expressions involving six covariant derivatives. It consists of all possible index 
contractions of 

1) [D^DMUDaDpD.DsU]., 

2) [D^DM]^ {[DaDfsD^UUDsU]- + rev) 6,^^, 

3) [D^DM]- {[D^DpUUD^DsU]^ + rev) e^,^,, 
4a) [D^DM]- i[DaDf3U]4D^U]4DsU]-+rev) , 
46) [D^DM]- i[D^U]4Df,D^UUDsU]^ + rev) , 

5) [D^DM]- {[D^UUDfsUUD^UUDsU]^ + rev) e,,,., 

6) [D^U]4D,U]-[DaUUD(,U]4D.,U]4DsU]-+Tev. (40) 

Here 'rev' stands for the terms in reverse order and it is understood that we must take all 
possible index contractions and in the end a trace. 

The general result for multiple traces will be essentially the same, except that we need 
to take all possible multiple traces, as well as all index contractions, and except that the 
reversed term in each of the above is modified so as to reverse the terms in each trace 
separately. 

We are now in a position to simply read off the results from Eq. (0), as the only thing 
that really needs to be done is to be sure that all independent index contractions are included. 
Note that a number of the possible contractions vanish because the multiple derivatives are 
symmetric in their indices while eajj^ys is antisymmetric. Also we are interested only in the 
structure, so overall factors will be neglected. 

For fx = u we get from 1) 

Tr [[D^.D'^UUD^D-DpD^U]^) , (41) 

from 4a) 

Tr {[D^D^UUD^D^^UUDpUUDfU]^) , (42) 
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Tr ([D^D'^C/]_([L»„D^C/]_[L>"C/]_[L»%]_ + [L>^C/]_[D-C/]_[D,L>^C/]_)) , (43) 
from 4b) 

Tr {[D^D^UUDpUUD^D^UUD^U]:) , (44) 

Tr {[D^D>'U]_[D''U]_[D^DpU]_[D^U]_) , (45) 

and from 5) 

Tr e^p^s. (46) 

Terms from 2) or 3) all vanish, using the symmetry of multiple covariant derivatives. 

For fx ly, recall that we cannot have repeated indices in the multiple derivatives. Thus 
we obtain from 4a) 

Tr{[D^D,UUD'^D''UUD^UUD"U]_) , (47) 
Tr {[D^DM]_[D'^D''UUD''UUD^U]_) , (48) 
Tr {[D^DMUD^D-UUD^UUD'-U].) , (49) 

from 4b) 

Tr {[D^DM]-[DaU]-[D'^D''UUD''U].) , (50) 

Tr {[D^D,UU[D''UUD''D''UUD^U]. + [D,t/]_[D^D"t/]_[D^C/]_)) , (51) 
and from 5) 

Tr {[D''D''UU[DM]-[D^UUDmUD'U]. + [D^t/]_[D^t/]_[D'3[/]_[D,t/]_)) e^p^s, 

(52) 

Tr ([T>'^T>^C/]_([T'^[/]_[T>,C/]_[T>^[/]_[T>^[/]_ + [T>^C/]_[T'^C/]_[T>,C/]_[T>^[/]_)) e^f^^g- (53) 
Finally from 6) we get 

Tr ([T>,C/]_[T>"C/]_[T>^C/]_[T'^[/]_[T'^C/]_[T>^C/]_) , (54) 

Tr {[D^UUD'^UUDf.UUD^UUDf'UUDm].) , (55) 
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Tr {[D^UUD'-UUDpUUD^UUD^UUD^Ul 



Tr ( [D„ U] _ [Dfs U] _ [D^ t/] _ [D" t/] _ [D^ t/] _ [Z^^ f/] . 



Tr ( [D„ U]4Df3U]4D^U]4D^U]4D''U]4Dm]- 



(56) 
(57) 
(58) 



A number of these terms differ from others only by a reordering of terms and thus we 
might expect that the trace relations of Sec. [11E6| might connect some of them. This is 
in fact the case. The sets Eqs. (||) and (H), Eqs. (||) and (||), Eqs. (0) and (0), Eqs. 

in each case contain all of the 



(H), (H), and (IT]), and Eqs. (H, (H), (H), (P, and (H, 
independent permutations. Hence one of each set can be eliminated using the trace relation, 
say Eqs. (0), (H), (0), (|5l]), and (H) 



Likewise Eqs. (|§), (|52D, and ( ^31) are related by 
the epsilon relation of Sec. |ilE7| and we can use it to eliminate say Eq. (0). 

It is interesting to note that this epsilon relation connects a term involving the equation 
of motion, Eq. (PB|), with two that do not obviously contain the equation of motion, Eqs. 
(^) and (|53|) . As we discuss later, the equation of motion terms can be transformed away 
by an appropriate transformation of the fields, and so do not contribute. Had we chosen to 
use the epsilon relation to eliminate the equation of motion term Eq. (|46|) rather than one 
of the others we would have ended up with two terms from this set instead of one, and it 
would have been extremely non obvious, from just looking at the terms, that one was in fact 
redundant. 

Consider now the multiple trace terms arising from the general result of Eq. (^). We 
simply have to construct all possible combinations of traces, plus all possible contractions of 
indices. The terms arising from charge conjugation have to be modified slightly, as we take 
the terms in reverse order within each separate trace. The result is simplified very much by 
the fact that Tr ([D^[/]_) = Tr {[D^D^U]-) = 0. That fact, plus the symmetry of multiple 
derivatives, means that there are no contributions from lines 1), 2), 3) or 5) of Eq. (|40|). 

With fi = ly we get from 4a) 



Tr ([D^D^'/7]_[D„D"[/]_) Tr {[D^UUD^Ul 



(59) 



and from 4b) 



Tr{[D^D''U]4D^D(,U]-)Tr([D''U]4D''U]- 



Tr{[D^D>'UUD^U]^)Tr([DpD''U]4D''U] 



(60) 



(61) 



Tr{[D^D^^UUD^U].)Tr[[D^D^UUDf,Ul 
With fi ^ ly we get from 4a) 

Tr{[D^DM]-[D^D''U].) Tr {[D^UUD^U]. 



(62) 



(63) 
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Tr{[D^DMUD''D''U]-)Tr{[D''U]4DaU]-) , (64) 

and from 4b) 

Tr{[D^DM]-[DaU]-)Tr{[D>'D''UUD"U]^) , (65) 

Tr{[D^DM]-[D''U]^)Tr{[D>'D''UUD^U]^) , (66) 

Tr{[D^D,UUDaU].)Tr{[D>^D^UUD''U].) . (67) 
Finally from 6) we get 

Tr{[D,UUD>'U].)Tr{[Do^UUD"UUD^UUD''U].), (68) 



Tr{[D^UUD^^U]^)Tr [[D^UUDpUUD-UUD^U]^) , (69) 

Tr{[D^UUDM]-)Tr{[D^^UUD''UUD^UUD"U]^), (70) 

Tr([D^[/]_[D.f/]_)Tr([D^f/]_[D,f/]_[OT]_[D-f/]_), (71) 

Tr{[D,UUD>'UUDo.U].)Tr([DpUUD''UUD"U].), (72) 

Tri[D^UUDM]-[DaU].)Tri[D'^UUD''UUD''U]^), (73) 

Tr ([D^f/]_ [D,t/]_)Tr([OT]_[D"f/]_[D''t/]_), (74) 

Tr {[D.UUD^^U]^) Tr ([D,f/]_ [D"f/]_) Tr ([D^f/]_ [D'^f/]^) , (75) 

Tr {[D^UUD^U]^) Tr {[D^UUD^U]^) Tr ([D"UUD''U]^ , (76) 



Tr {[D^UUDM]-) Tr {[D^UUD^U]^) Tr {[D'^UUD^U]^) . (77) 

There are two trace relations for the multiple trace terms, one relating Eqs. (|68|) and 
3) and one relating Eqs. (|70|) and (ffll). We use these to eliminate Eqs. (|69| ) and ([711). 
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B. Terms with four L>^'s and one x'^'^ 

Terms at this level will have four D^s and one x^'^ 7 which may be a simple x if /i = 
or may be G'^'^ or H^'^ constructed from the F'^'^^s according to Eq. (pOD if fi u. In the 
former case there will be four Lorentz indices which may be contracted in two pairs or with 
an eafB-yS symbol. In the second case there will be six indices and the types of contractions 
will be the same as in the preceding section, except that x'^'^ is antisymmetric in its indices. 

Proceeding as before we write down schematically the types of terms which are possible, 
using the initial building blocks, with the understanding that eventually traces of these 
terms will be taken. We find 

1) {D'U)iD\n, 

2) {D^U){D^U)x'"'U, 

3) {D'U){DU){DU)x''' , 

4) {D^U){DU)[DxnU, 

5) {D\^''){DU){DU)U, 

6) {Dxn{DU){DU){DU), 

7) x^''U{DU){DU){DU){DU). (78) 



In writing these down we have already used the total derivative arguments of Sec. [II E 1 
to eliminate terms involving D'^ and D^. In fact it is a general result that one needs at most 
half of the covariant derivatives in the term acting on a single factor. We have included 
explicitly the minimum number of t/'s necessary to give an even number of factors so that it 
is possible to get the right overall transformation properties. Note also that at the present 
stage there may in principle be additional nonadjacent U and W pairs distributed through 
some of the terms. Using again the total derivative argument one of the D's in 5) can be 
moved off of the x to give terms like 6) and 4). Similarly in 6) the D can be moved from 
the X to give terms like 3) and, if there are extra U,W pairs, like 7). Thus 5) and 6) need 
not be considered. 



Using Sec. [11E2| all multiple derivatives can be taken to be symmetric in their indices 



and this allows us to simplify 1) and 2). In these terms, to preserve Lorentz invariance, at 
least two D's must be contracted together, or contracted with an eap'ys- Using the index 
exchange result of Sec. [1IE3| these two D's can be brought together acting on a U, with 
the additional terms produced being either total derivatives or like some of the other terms. 
By symmetry the eap-ys term must vanish. Likewise by symmetry, the remaining two D^s 
cannot be contracted with the x'^'^ ^ so must also be contracted together. Thus x^'^ must 
be simply x aiid we find for 1) and 2) 

1) ^ {D^D''U){Df,Df'x), 

2) ^ {D^D^U){DpDPU)xU. (79) 

As before to get an expression transforming as Vr . . . V^j, so that the trace will be invari- 
ant, we write each factor in Eq. (^) in the form [A\±. Then Sec. |II E 4| and particularly Eq. 



(p9|) allows us to eliminate [D [/]_|_ in favor of terms we already have. 
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The result of all these manipulations is the set of general forms 

1) [D^D'^UUD^D^xh, 

2) [D^D'^UUDpD^UUxh. 

3) [D^DpUUD.UUDsUUx'''']^, 

4) [D^DpUUD,UUDsxn±, 

7) [xn±[D^UUDpU]4D,UUDsU]^, (80) 

where it is understood that we still must take all possible contractions of indices, perhaps 
with eQ,/3-y5, and all different independent orderings of terms. 

The last step is to impose parity and charge conjugation invariance. In accord with Sec. 
II E 5| requiring parity invariance simply leads to the overall factor 1 + (—1)*+^+'= just as 
before. P is the sum of intrinsic parities p, but in effect here just counts the number of 
H'^'^^s, as all other quantities have p = 0. For the five terms in Eq. ([80|) these factors are 



(1 T 1), (1 ± 1), (1 T (-1)^+'), (1 ± (-1)^^'), and (1 ± (-1)^+'), respectively, where the 
upper (lower) signs correspond just to the factors [A]^ ([^]-) involving x or x^'^ ■ Similarly 
to get charge conjugation invariance we must add a term (— l)'^x the original term with all 
the major factors in reverse order. 

We can now put all of this together to get the general result for the four D one x^" 
terms analogous to the six D result of Eq. (|iO|). It is most convenient to separate the 
/J, = u, i.e. X, terms from the 7^ z/ terms as there are some 'accidental' simplifications 
which are different in the two cases, and since the values of the intrinsic parity and charge 
conjugation quantum numbers are different. For ^ = u observe that the symmetry of the 
multiple derivatives ensures that there will be no ea/s-ys terms arising from 3) or 4) in Eq. 
(^). For jj, ^ u we have by explicit calculation [x'^'']- = 0. Furthermore we can show, by 
taking the covariant derivative of [x'"^]- and using the chain rule Eq. (|T^ that [D^x'^'^J-can 
be expressed in terms of other factors we have kept. Hence it is not independent and can 
be dropped. 

We thus find the most general form for single trace expressions with four D's and one 
X^'^ to be, for fx = u 

1) [D^D^UUD^D^x]^, 

2) [D^D-UUD^D^UUxU + Tev, 

3) [D^DpUUD,UUDsUUx]-+rey, 

4) [D^DpUUD.,UUDsx]+ + rev, 

7a) [xUlDaUUD^sUUDyUUDsU]- + rev, 

7b) [xUDaUUDfsUUDyUUDsUW^^' + rev, (81) 



and ioT fi ^ u 



3a) [D^Df,UUD^UUDsUUG''''Ueprr,a - rev, 

3b) [D^Df3U]4D^UUDsUUH'''']+ + rev, 

4a) [D^Df.UUD^UUDsG^^'U - rev, 

46) [D^DpUUD,UUDsH>^''Ueprn„ + rev. 
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7a) [G"^1 + [D«f/]_[D^[/]_[D^?7]_[D5[/]_ - rev, 

7b) [Hn4D^UUD^gUUD,UUDsUUpr,^ + Tev. (82) 

In both of these results 'rev' stands for a piece with the factors in reverse order and it is 
understood that we still must take all possible contractions of free indices and all possible 
independent orderings of factors and, to get the final contribution to the Lagrangian, take 
the trace. 

Everything done to derive these results holds also essentially unchanged for multiple 
traces, by virtue of the total derivative arguments of Eqs. (P^ ) and (^). We simply must 
take all possible combinations of traces, instead of a single trace, and must reverse the order 
of the terms in each individual trace separately. 

It is now straightforward to tabulate the terms arising from Eqs. (|8T|) and (|82D . Starting 
with Eq. (^) we get for the terms involving x from 1) 



from 2) 



from 3) 



Tr([D„D°[/]_[D^D^x]-), (83) 



Tr ([D^D^UUDpD^UUxU) , (84) 



from 4) 



Tr {[D^D-UUlDpUUD^UUx]- + [x]-[DpUUD^U]^)) , (85) 

Tr{[D^D^UUDfsUUx]-[D''U]^) , (86) 

Tr {[D^D,UU[D-UUD^UUx]- + [xUD^UUD^U]^)) , (87) 

Tr([D,D^[/]_[D"[/]_[x]-[D^[/]-) , (88) 

Tr {[D^D^UU[Df,UUD^x]+ + [D^xUD^U]^)) , (89) 

Tr ([D^DfsUU[D^UUD^x]+ + [D^xUD'^U]^)) , (90) 



and from 7) 



Tr {[xU[D..UUD"UUD^UUD^U]^) , (91) 
Tr {[xUDaUUD^UUD^UUD^U].) , (92) 
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Tr ([xUDaUUDpUUD^UUD'^U]^) , (93) 



Tr{[x]-[DaUUDpUUD,UUDsU]_)e"^'''. (94) 

It is equally straightforward to read off the multiple trace terms originating from Eq. 
Bl|). Recall that Tr {[DaDpU]^) = Tr ([Dq[/]_) = 0. Thus we have no terms from 1) and 

Tr{[D^D^UUDpD^U]^)Tr{[xU), (95) 



get from 2) 



from 3) 



from 4) 



and from 7) 



Tr{[D^D^UUDf,U]^)Tr[[D''UUx]-) , (96) 

Tr{[D^D^UUx]-)Tr{[DpUUDfU].) , (97) 

Tr {[D^D^UUDpUUD^U]^) Tr {[x]-) , (98) 

Tr {[D^Df,UUD-U].)Tr {[D^UUx]-) , (99) 

Tr {[D^DpUUx]-)Tr ([D-UUD^U]^) , (100) 

Tr ([D^DpUUD'^UUDf'U]^) Tr {[x]-) , (101) 

Tr{[D^D'^UUDpU]^)Tr{[Df'xU) , (102) 

Tr{[D^DpUUD'^U]^)Tr{[D^xU) , (103) 

Tr {[x]+)Tr ([D^UUD^UUDpUUD^U].) , (104) 

Tri[xh)Tr([D^UUDpUUD'^UUD^U]) , (105) 



Tr{[x]4DaU].)Tr{[D-UUDpUUD^U].) , (106) 
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Tr{[x]4DaUUD''U]^)Tr([DpUUD^U]^) , (107) 
Tr {[x]4DaUUDpU]4Tr [[D^UUD^U]^) , (108) 
Tr([x]+)rr([D,[/]_[D"[/]_)rr , (109) 



Tr{[xU)Tr{[D^UUDpU]4Tr[[D^UUD^U]^) . (110) 
Like in the previous section there are here also several sets of terms consisting of all 



permutations of the factors. Thus we can use the trace relations of Sec. |II E 6| to express 
one of each set in terms of the others and of terms with more traces which are lower down 
in the hierarchy. In particular Eqs. ( pSj) and (|86D, Eqs. (|^) and (|5), Eqs. (|9T|), (0), and 
(P3|) and Eqs. ( |104D and ( |105D form such sets and we choose to eliminate Eqs. (P^D, (P7D, 
(0), and ( |105|) using the trace relations. 

We now evaluate the terms coming from Eq. (|82|) which contain G^'^ and H^'^. We obtain 
for the single trace terms from 3a) 



Tr [[D^D"UU[D''UUD^UUGn+ " [Gn + [D"'UUD''U]4) e^,,,, (111) 

Tr {[D^D-UUDf'UUGn+iD^U].) e^^^,, (112) 

Tr ([Z}"D^^]_([Z}^f/]_[D^f/]_[G'n+ - [G'^'UD^'HUD^Uj^) e„^^,, (113) 

Tr {[D'^D^UUlDpUUGn + lD^'U]- - [DmUGn + lDpU]^) e„^^,, (114) 

Tr {[D'^D^UU[D''UUDf,UUGn+ " [G^'l + i^/^^]- [^"f^]-)) W> (US) 

Tr ([Z}-D^t/]_([D^t/]_[D^[/]_[G"^n+ - [G^'UD'UUDm^)) e„^,,, (116) 

Tr ([D"D^t/]_[D^t/]_[G'^1 + [D^f/]_) e„^,,, (117) 

from 3b) 

Tr ([D,D"t/]_([D^t/]_[D,t/]_[//'^-]+ + [H'^'']4D,UUD,U]4) , (118) 



Tr [[D^Df,UU[D^UUDM]-[Hn+ + [i^H + f^-^^]- [^^^f^]-)] , (US) 
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Tt [[D^DpU]_{[DM]-[D^UUHn^ + + , (120) 

Tr {[D,Df,U].{\D,UUH>^']+[D^U]. + [D^UUH^^] + [DM]-)) , (121) 

from 4a) 

Tr {[D^D^U].{\D^UUD,G^-]+ - [D,G^^]+[D^U].)) , (122) 

Tr ([D.D^UUiD^UUD^n^ - [D^Cn + iDpU].)) , (123) 

Tr {[D^D^UU[DM]-[D/3G''''U - [^/?G"^1+[^.f^]-)) , (124) 

from 4b) 

Tr {[D^D^UU[D^UUD^H'^''U + [D^H^^D^UU)) ep,,,, (125) 

Tr ([D"D^t/]_([D^t/]_[D^i/n+ + [D^H^^'UD^U].)) e^,,,, (126) 

Tr {[D^D^UU[D'^UUDpH^^^]+ + [DpH^^]^[D'^UU) e^,,,, (127) 

Tr [[D''DpU].{[D^U].[D'H^^]+ + [D' H^-] + [D^U].)) e«^5., (128) 

from 7a) 

Tr {[G^'']+{[D,U]_[D,U]_[D^UUD-U]_ - [D^U]_[D'^UUD,UUD^U].)) , (129) 

Tr {[G^']+{[D^U].[D^UUD,UUD'^U]. - [D^U].[D,U].[D'^UUD^U].)) , (130) 

Tr {[G^'] + [D^U].[D^U].[D^U].[D,U].) , (131) 

Tr {[Gn + [DaU]-[D>.UUDM]-[D"U]^) , (132) 

and from 7b) 

Tr {[H^^'MD^UUDf'UUD.UUDm]. + [D^t/]_p^t/]_p'3f/]_[i?"[/]_)) e^.«^, (133) 

Tr {[H^^'^MD'^UUD^UUD^UUD'^UU + [D^C/]-[L'^C/]_[L>^[/] _[£>"[/]_)) e^^^^, (134) 
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Tr [[m''] + {\DM]_[D''UUD'<U]4D'U]^ + [D^f/]_ [D^?7]_ e,p,s.. (135) 



Tr[\W^''] + {\DPU].[D,UUD'<UUD'U]^ + [D'UUD'<UUD,UUD('UU) t^p^s- (136) 



In this group the set Eqs. ( tL19| ), ( |120| ), and (|121|) are related by the trace relations of Sec. 
ilE6| , and we use that result to eliminate Eq. (|121| ). Also a number of terms are related by 



the epsilon relations of Sec. |II E 7| . In particular there is one relation among terms of each 
of the sets Eqs. (|TT2|), (|n|), and ([n^), Eqs. (|m]), (^j, (|Tl5]), and (^) and Eqs. ( pSj ), 
( |126|) , ([12^), and (|128D which we use to eliminate say Eqs. ([TT|), (|l|), and (112^). Just as in 
the discussion following Eq. ( pBD we here also keep the equation of motion terms, since they 
can be eliminated later via field transformations, and use these relations to eliminate other 
terms. There are two relations among the set Eqs. ( 133 ), ( 134 ), ( 135 ), and (136) which we 
use to eliminate Eqs. ( |135| ) and (|136|) . 

Again we can write down directly from Eq. (|8^) the multiple trace terms. Many possible 



terms vanish because for /i 7^ z/ we have Tr ([x 



flUl 



0, as can be shown 



by explicit evaluation using the definitions Eq. (|20|) . Because of these relations, and also 
Tr {[DaDpU]^) = Tr {[DaU]^) = 0, we obtain no terms from 3a), 4a), or 4b), and we get 
from 3b) 



Tr {[D^D'^UUD^U]-) Tr {[DM]-[H' 



IJ.U] 



(137) 



Tr([D^DpUUD''U])Tri[DM]- [H 



Tr {[D^DpUUDM]-) Tr [//''^] 



J+ ' 



Tr {[D^DpUUH^^-'U) Tr ([D^UUDM] 



from 7a) 

Tr {[Gn+[DaU].)Tr {[D'^UUD^UUDM]-) , 
Tr{[G^'']4D,UUDM]-)Tr{[D^UUD^U]_), 

Tr([Gn + ([I>M^]-[^a^]- - [D^UUD^U].))Tri[D''UUDMl 
and from 7b) 

Tr{[H'^''UD''U]^)Tr([D^UUD.,UUD'^U])e^,^p, 



(138) 
(139) 
(140) 

(141) 
(142) 
(143) 

(144) 



Tr{[H^''U{[D''UUDm]^ + [DmUD'^U]^))Tr{[D^UUD,U]^)e,,^p. (145) 
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C. Terms with two L>^'s and two x^'^'s 

Consider now terms of the third level in the hierarchy, those containing two -D^'s and two 
X^'^'s. The x'^^'s can be just x if = or can be the full involving the F'^'^^s if ^ u, 
or there can be one of each. Because of these different possibilities, and the fact that each 
X'^" has its own parity and charge conjugation quantum numbers there are so many options 
that it is much more difficult in this case to develop a general result and simply evaluate 
it, as was done in the previous two sections. We will however proceed as far as possible in 
general, but will be forced in the end to simply enumerate the possibilities. 

Note that we will have to distinguish between the two x'^'^'s. To do this we will call one 
X and the other x- To simplify the notation we will also drop the Lorentz indices on x until 
the end. 

We can then write the possible terms directly in terms of the building blocks [A]± as 



1) 


[DX]± 






2a) 


[DUl 




[X]: 


2b) 


[DU]- 


mu 


M: 


3) 


[DU]- 


[DU]- 


Ix] 



X]±- (146) 



In these expressions we must take all independent orders, and of course eventually put the 
Lorentz indices back in and take all possible contractions. The ± signs on the x and x terms 
are uncorrelated, so that all combinations must be included. 

To obtain these expressions we have already systematically used the total derivative 
argument of Sec. |11 E 1| to remove any D^'s. However this procedure may have hidden some 
equation of motion terms in this set, which we want to extract. For example, moving the 
D acting on x in 2a) over to the DU using the total derivative argument could generate a 
[DaD°'U]- term. There are some subtleties in such a transformation however. So far we 
have eliminated terms high in the hierarchy in favor of those lower down, i.e. in this case 
we eliminate those with more D's acting on a single factor in favor of those with fewer D's. 
Since we always keep the most general set of terms at the lower level, it is never necessary 
to actually work out explicitly the relation used to eliminate the higher level term. Here 
however we are asking whether a specific lower level term can be eliminated in favor of a 
higher level one. Thus to be sure one is not over or under counting the independent terms 
it will be necessary here to work out the transformation explicitly. 

To do this consider the total derivative {[DuU]_[x]±[x]±) , where the argument of the 
derivative is understood to be a single or multiple trace. Since total derivatives can be 
dropped in the Lagrangian, this can be treated as if it were effectively zero. Thus using the 



results of Sec. [II E 1| and particularly Eq. (|2^) we find 

= [D,DM]-[X]±[X]± 

+ [DMUD,X]±[X]± + [DM]-[X]±[D,X]± 

+ terms of the class 3), (147) 

where again the traces have not been put in explicitly. The terms of the class 3) arise from 
the covariant derivative acting on the U and W in the various factors [A]± together with 
some algebraic rearrangement to express everything in terms of the standard factors [A]±. 
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Thus we see that we can express a particular sum of the terms from the classes 2a) and 
2b) in terms of a class containing a double derivative term [D^D^U]^. To get this particular 
sum we must take a new basis for the terms of 2a) and 2b) consisting of the sums and 
differences of terms. When this is done we can use this double derivative term instead of 
the sum if we wish, but we must keep the remaining (difference) terms from 2a) and 2b). 
Note that we could have started this discussion with alternative ordering, x but since 
we always include all orders of the factors, this gives the same results. 

This means that in the general result of Eq. ( |146|) we may use instead of 2a) and 2b) 

2a)' [DU]^i[Dx]±[x]±-[x]±[Dx]±), 

2by [D^DMUxUlxU, (148) 

where as usual, all orderings must be taken and the ± signs are not correlated. 

Observe that everything done so far applies to multiple traces just as to single traces. 
One sees, by following through the arguments leading to Eq. (148), that if we start in 2a) 



and 2b) with say, the trace of the first two factors times the trace of the third, then we will 
also get the trace of the first two factors times the trace of the third in each of the terms 
making up 2a)' and 2b)'. 

Finally it is possible to apply essentially the same argument to 1) of Eq. ( |146[ ) and show 
that we can use either of the alternative forms 

1)' [xUD^D^xh, 

1)" [D,D,x]±[x]±- (149) 

Next we must ensure that parity and charge conjugation are satisfied. Recall that effec- 
tively under parity [A]± —>■ ±{—iyW[A]±U with p being the intrinsic parity of the operator 
A. There is also the extra factor of (—1)*^, where e counts the number of e^/j^^'s in the term. 
Let p and p be the intrinsic parities of x and x let (—1)^ here account for the product 
of ± signs coming from the [x]± and [x]± factors only, e. g. s = for + and s = 1 
for Thus we get a parity invariant term by multiplying the forms 1) and 3) by 

(1 + l-iy+P+P+') and the forms 2a)' and 2b)' by (1 - (-1)«+p+p+^). 

To ensure charge conjugation invariance we must add (—1)^''"'^ times the major factors 
in reverse order, where c and c are the intrinsic charge conjugation quantum numbers of 
X and X- For terms of type 1) and the multiple trace terms from type 2a)' and 2b)' with 
only two factors in the trace the cyclic property of the trace means that this amounts to 
multiplying by the overall factor (1 + (—1)'^"'"'^). For the others the reversed term has to be 
added explicitly, though in many cases adding such a term makes different initial starting 
orders for the factors give the same result. 

We can now summarize the general form for the two and two x^'^ terms as follows. 

1) (1 + (-i)^+^'+^^+^)(i + {-ir')[D,x''%[Dur']±, 

2a)' (1 - (-l)WP+e) [[D,UU[D,x''^Ux'']± - [x''^]±[D.r']±) + (-1)'=+' (rev)) , 
26)' (1 - (-l)WP+e) [[D,D,UUx''^Ur']± + (-1)^^' (rev)) , 
3) (1 + (-1)^+^^+^+^) [[D,UUD,UUx-^]^[x'']± + (-1)'=+' (rev)) . (150) 
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In these expressions one must contract the Lorentz indices in ah possible ways, including 
perhaps contracting with an e^g^^, must take all possible independent orders for the factors, 
and in the end must take a trace. The ± signs are not correlated, though the (—1)* in the 
parity factor will enforce a correlation for given values of the other quantum numbers. The 
indication 'rev' means to take the major factors in reverse order. For multiple traces, one 
must in addition take all possible combinations of different traces. Just as in the previous 
section the evaluation of the result, particularly the multiple trace part, is simplified by 
the relations Tr {[D^D^U]_) = Tr([D^[/]_) = 0. Also for 7^ Tr([x^'']±) = [x'^l- = 
Tr {[DaX^'']+) — 0. Furthermore [DaX^"]- can be expressed in terms of other quantities 
lower in the hierarchy which have been kept, and so can be dropped. 

We now proceed to evaluate this general expression. There are three main cases, namely, 
I) X s-iicl X both simple x'Si H) one a x and the other a x^'^ with fi ^ u, and III) both 
^/ii/5g Furthermore each x'^" can be a G'^'^ or an H'^". We will evaluate the single and 
multiple traces together, as that seems simplest here, and will always drop irrelevant overall 
numerical factors. 

Consider case I) where both x and x ai'c simple x's. The quantum numbers p = p = c = 
c = e = so that s — for terms from 1) and 3) and s — 1 for those from 2a)' and 2b)'. We 
obtain from 1) 

Tr{[DMD^xU): (151) 



Tr{[D,xUD'x]-), 



Tr{[D^xU)Tr{[D^xU): 



Tr{[D,x]-)Tr{[D^x]-), 



from 2a)' 



Tr{[D,UU[D^x]4x]- - [xUiD'x]- + [xUD'xU - [D'xUxU)) , 



Tr {[D.UUD^xU) Tr ([x]-) - Tr {[D.UUxU) Tr {[D^x]-) , 



from 2b)' 



Tr {[D^UUD'^x]-) Tr ([x]+) - Tr {[D.UUx]-) Tr {[D'^x]+) , 



Tr{[D,D'^UU[x]4x]- + [xUxU)), 



Tr{[D^D^UUxU)Tr{[x]-): 



Tr{[D^D^UUx]-)Tr{[x]+), 



(152) 
(153) 
(154) 

(155) 
(156) 
(157) 

(158) 
(159) 

(160) 



29 



and from 3) 



Tr{[D,UUD^UUx] + [x]+), 



TrilD.UUxUD'-UUxh), 



Tr{[D,UUD>^UUx]-[x]-), 



TrilD^UUxUD'-UUx]-), 



Tr{[D^UUD'^UUx]+)Tr{[xl 



Tr{[D,UUD^UUx]-)Tr{[x]- 



Tr{[D,UUD>^U].)Tr{[xU[x]- 



Tri[D^UUD'^U]^)Tri[xUxl 



Tr{[D^UUx]+)Tri[D'^UUxl 



Tri[D^UUx]-)Tri[D'^UUx]- 



(161) 
(162) 
(163) 
(164) 
(165) 
(166) 
(167) 
(168) 
(169) 
(170) 



Tr {[D.UUD^'UU Tr (M+) Tr (M+) , 



(171) 



Tr ( [D^U] _ [D'^U] ^)Tr{[x]-)Tri[x]-). 



(172) 



Of these there are two sets of terms which consist of all permutations of factors, Eqs. ( |161| ) 
and (|162| ) and Eqs. ( p.63| ) and (|164|) . Within each set the terms are thus related by the trace 
relations of Sec. |11 E 6| which we use to eliminate Eqs. ( |162D and ( |164| ). 



Now consider case II) where say x is the simple x ^md x is one of G'^'^ or if^*^. Now 
p = c = 0, with p = and c = 1 for G^^ or p = 1 and c = for H'^^. There are now four 
indices so it is possible to have an efj^^^ys factor. Observe that 2b)', and also 1), if we use the 
equivalent form 1)' of Eq. ( [L49| ), are symmetric in the interchange fi u and hence they 
will vanish when contracted either with x^'^ or e^i^^s- Thus non zero terms come only from 
2a)' and 3). 

We then obtain from 2a)' without an e^^^s factor 



Tri[D,UU[D,x]-[Gn+ " [x]-[D.G^'U " [G^AD^x]- + [D^GnAx]-)) , (173) 
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Tr {[D^UUiD^xUH'^U - [xUD.H^^U + - [D^H>^^UxU)) , (174) 

Tt {[D.UUm^]^) Tr ([D.x]+) - Tr {[D.UUD^H^^]^) Tr ([x]+) , (175) 
and with an e^^^s factor 

Tr [[D^UUlD^xUH-^'U - [xUD^'H^'U + [H^\[D'x]- - [D^H-^'Ux]-)) W> (177) 

(Tt {[D>^U]_[H^\) Tr (p^x]-) - Tr {[D>^U]_[D^ H^\) Tr ([x]-)) t.^^s- (178) 
Prom 3) without an e^„rys factor we get 

Tr{\D,UUDMU[x]AGn+ + [G"^1+[x]+)) , (179) 

Tr([Z),t/]_[x] + [D.t/]_[G'n+), (180) 

Tr ([D,[/]_[D.t/]_[G"^'^]+) Tr ([x]+) , (181) 

Tr ([D,t/]_[D.t/]_([x]_[i/n+ - [HnAx]-)) , (182) 

Tr Tr {[DM]-[x]-) , (183) 

and with an e^^-ys factor 

Tr {[D^^UUD^UUixUG-^'U + [C^'ln-M-)) W> (184) 

Tr ([D'^t/]_[x]-P'^t/]-[G'^']+) (185) 

Tr ([D^UUD'^UUG^'U) Tr ([x]-) e^.^^, (186) 

Tr ([T>''C/]_[i;'^C/]_([x]+[//^^]+ - [iy^1+[x]+)) e^.,., (187) 

Tr Tr ([D^[7]_[x] + ) e,,,s. (188) 



Consider now the third case in which both x and x are x'^'^ with /j, ^ u. Observe first 
that s = since aU factors involving x^'^'s are of the form [. . •x'*'^ • • •]+• This means that 
when p + p + e is even only 1) and 3) contribute while when it is odd only 2a)' and 2b)' 
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contribute. Furthermore for the HH or GG terms p + p + e e and c + c is even whereas for 
the GH terms p + p + e^l + e and c + c is odd. All of the individual factors are traceless, 
so the only multiple traces possible must involve two factors and thus can come only from 
3) . Finally by considering various permutations we see that we get fewer terms at this level 
if we use the 1)' form instead of the 1) form. 
We thus obtain from 1)' 




(189) 




(190) 




(191) 




(192) 



From 3) we get the following terms without an Cap^ys 




(193) 




(194) 




(195) 



Tt [[D^UU[G-'^UDPuUGp,]+ + [Gp^UD^UUG^^U)) , 



(196) 




(197) 




(198) 




(199) 




(200) 




(201) 
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Tr([D,UUD^UUH^%[H^pU,, 



/3J+ ' 



(202) 
(203) 
(204) 



Tr ([D^UUD^UUH^-'UHp,] 



Tr {[D^UUD^UUHMH'-^U) ^ 



Tr{[D^UUD^U].)Tr([H^P]4H^ 



(205) 
(206) 
(207) 
(208) 



Tr[[D,UUH^PU)Tr{[D^UUHo.p\ 



\+) 1 



Tr ([D^UUD^U]^) Tr {[H^^UH^.U) , 



Tr {[D^UUH"''U) Tr ([D^f/]_ , 



Tr {[D'^UUH^^'U] Tr{[DpUUH^,] 



(209) 
(210) 
(211) 
(212) 



Of these there are four sets of terms whose elements are related by the trace relations of Sec. 
PE^ , namely, Eqs. (|193]) and (|9|), Eqs. (|9|), (|196D, and (|193), Eqs. (|03D and (|20|), and 
Eqs. (|205|), (|206D, and (|207|). We use those relations to eliminate one term of each set, Eqs. 



O, dH), dH, and (11). 

The eais-ys terms from 3) are 



Tr ([D.UUD'^UUiG'^f'UlH^'U " [i/^'] + [G°^]+)) e„ 



1/375 ; 



(213) 



71 r/C"/3l 



(214) 
(215) 
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Tr 












(216) 


Tr 








-[H^']4D^UUG^%] 




(217) 


Tr 








- [D^UUH^'UG^^U] 




(218) 


Tr 




-{[DaU] 




- [g-^'uIh^Wd^uu 




(219) 


Tr 








- [G^'^UiD^UUH^^U] 


j ^M/375' 


(220) 


Tr 


{[D^U]^ 






- [D^UUG^'UIH'^^]^ 




(221) 



There are two sets of epsilon relations, as described in Sec. PE7| , among this group of 
terms. One consists of two relations among the set Eqs. ( |^15| ), ( |217| ), and ( ^2CI| ) and will be 
used to eliminate Eqs. (|217|) and (|220| ). The other involves three relations among the set 
Eqs. (|13D, (|1|), (^, (|ll), dTI), and (|2lD and will be used to eliminate Eqs. (gTsD , 



m, and M. 



From 2a)' we obtain 



(222) 



(223) 



Tr {[Dm]^{[D^G"^]4H,f,U " [G"^]ADaH,f,U " [H,f,U[D^G''% + [D^H,f,]4G 



(224) 
(225) 



(226) 



(227) 
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Tr ([D^/7]_([D'^G°^'^] + [G',^]+ - [G'"^] + [D'^G„^]+)) e^,^^, (228) 



(229) 



/3l 



(230) 



There are two sets in this group, Eqs. ( ^25]) , (|2|), and (^ and Eqs. (|228|) , (p29D , and 



30|) , each related by a single epsilon relation which we use to eliminate Eqs. (|227|) and 



Finally from 2b)' we get 



Tr {[D.D'^UUIG^Wh^pU - [H^pUG'-^U)) , (231) 

Tr {[D^D^UU[G^^]4H,^U " [^7^«] + [G'"1+)) , (232) 

Tr {[D^D'^UUH^^]4H^''U) e^^.s, (233) 

Tr {[D^D^UU[H"^]4H'''U + [H''']4H^^U)) e,p,s, (234) 

Tr ([D^D''^]_[G"^] + [G'^^]+) e^^^^, (235) 



Tr ([D'^D„[/]_([G"'3]+[G^^]+ + [G'^^]+[G'"^]+)) e^^^^,. (236) 



There are two epsilon relations for this group, which allow us to eliminate Eq. (|234|) in favor 



of Eq. ( |233D and Eq. ( |236[ ) in favor of Eq. ( ^35| ). Note that this is again a case where it 
is important, if we are to minimize the number of terms, to use the relations to eliminate 
the non equation of motion terms, as the equation of motion terms can be eliminated in a 
different way. 

D. Terms with no D^^s and three x'^'^'s 

Consider now the final case with three x'^'^'s and no covariant derivatives. This is rela- 
tively simple compared with the previous cases and we can write the general case directly, 

(1 + i-iy^^^^^^^^^^^)[x't''U[xf]±[xI% + i-ir~'''-''nxl%[xf]±)- (237) 
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Here pi and Cj are the intrinsic parity and charge conjugation quantum numbers of the three 
x's. The ± signs are not correlated and as before (—1)'* is the product of the signs coming 
from the parity transformation on the individual [^]+'s, counting + for [A]^ and — for 
[A]_. By virtue of the second term with factors in reverse order, which comes from charge 
conjugation, the interchange X2 and xs gives, up to a sign, the original expression. Thus 
the three x's can be treated as distinguishable from the beginning. 

For the simplest case with three simple x's all of the pi and Cj as well as e are zero, which 
requires s = also. Thus we get 

Tr{[x]MM+), (238) 



Tr{[xUx]-[xU), (239) 

Tri[xU)Tri[xh[xU), (240) 

Tri[xU)Tri[x]-[x]-), (241) 

Tri[x]-)Tri[x]-[x]+), (242) 

Tri[xU)Tri[xU)Tri[xU), (243) 

Tr{[xU)Tr{[x]-)Tri[x]-). (244) 

There are no terms with two x's and one x^'^ as there is nothing with which to contract 
the indices. With only one simple x the general form of Eq. ( ^37]) reduces to 

(1 + i-iy^^^^^^^nixMxfUxl'U + i-ir-'-nx'/UxfU) (245) 

since [x'^'^]- = for ^ u. We see from this that if X2 = Xs, so that both are if^'^'s or 
both are G'^'^'s, then s = 1 if there is an epsilon term and s = if not. This is reversed if 
X2 7^ Xs- The number of multiple traces will be quite limited because Tr {[x^'^]±) = when 
fj, I'. Thus we obtain 

TrilxUHnAHM, (246) 



Tr{[x]4Gn4GM, (247) 
Tr{[xU[Hn4G,uU~[GMHn+)), (248) 
Tri[xU)Tri[Hn + [HM, (249) 
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TT{\x]+)Tr{\Gn4G,A+). (250) 
Tr[[x]-[H''^UHA+)e.MS, (251) 



Tr(M-[G'"^] + [G^^]+)e„^,,, (252) 

Tr (M + ([^^°^'] + [G^']+ - [G^ + W'^'W) e.p,s, (253) 

Tt{[x]-) Tr + (254) 

Tr ( M - ) Tr ( + [G^''] +) e^^^,. • (255) 



Finally when all three x's are x^^^^ the number of terms is very limited because of the 
parity and charge conjugation factors and because all factors are traceless. We find 

Tr([Gn + [G^J + [G.1+), (256) 
Tr([G^-] + [i/J + [i/,"]+), (257) 

Tr ([Gn + ([G"'^] + [^/]+ - [H-^UGp-'U)) e,..,. (258) 



The last of these, Eq. ( 258 ), is identically zero by virtue of the epsilon relations of Sec. IIE7 . 



IV. SIMPLIFICATION AND REORGANIZATION OF TERMS IN THE 
LAGRANGIAN - FINAL RESULTS 



We have now derived in Sees. [Ill A| - [III D| the complete set of terms contributing to the 
order Lagrangian Cq. The results are scattered through these sections in the order they 
were derived. We now want to collect those results in one place in a form of a Lagrangian 
with effective coefficients analogous to the standard Gasser-Leutwyler Lagrangian. In the 
course of doing this we want to simplify the forms as much as possible and to reorganize 
them so as to select out those terms which are likely to be most immediately useful. 

A. Equation of motion terms 

In the course of the derivation we have extracted as many terms as possible which are 
proportional to the factor [D ^D^U]^ which we have called the 'equation of motion terms' 
and have stated, but not proved, that these can be transformed away. Details of this 
transformation procedure are given with respect to Cq in HS] and with respect to lower 



37 



orders in a number of earlier papers [0,^,0. Here we just outline the general idea, as it 
produces a large reduction in the final number of terms which need to be considered for the 
general Lagrangian. 

The lowest order Lagrangian is given by 



(259) 



From this one can obtain the lowest order or classical equation of motion O^qj^^ = where 

1(2) _ orn n^rn or.,i 1 \ (260) 



O'Zm = '2[D,D^U]^ - 2[x]^ + -Tr {[x]-) . 



^EOM 6ach term 



We observe first that we can make the replacement [D^D'^?7]_ 

where it appears, where as always we have dropped the irrelevant numerical factor. Since 

(2) 

we have the most general form the extra terms added and subtracted to get O^q^j are just 
terms we already have. 

Now make a transformation on the fields of the form 



U ^U' = exp{iS)U 



(261) 



where S = S"^ and Tr (S) = 0. It is a general result |44-^ that such a transformation 
does not affect measurable quantities such as the S-matrix. Applied to £2 it generates a 



correction to lowest order in S of the form 



SC. 



P2 

^0 



Tr (iSO 



(2) 

EOM 



(262) 



and we can choose an S of order so that this term cancels the equation of motion terms in 
£4. This transformation generates corrections at order as well both from the second 
order in S correction to £2 and the first order in S correction to £4. For our purposes these 
corrections can simply be absorbed in the terms of Cq since we have the most general form. 
Finally we make a second transformation on £2 using an 5* of order and thus generate 
a correction term analogous to that of Eq. ( p62| ) which is of order and proportional to 



O 



(2) 

EOM- 



(2) 

By choosing S properly we can eliminate those terms of Cq which contain O^Qjy,j. 

The conclusion one draws from this discussion is that in order to generate the most 

(2) 

general Cq in its simplest form, we can just drop all terms proportional to O^q^j. This 
allows us to ehminate the 23 terms given by Eqs. (|T]), (g^), (||), (§Up, (^, (^), (§2|), 

(H), (n, (p, dn), m, ®, (0), (0), o, dni), (imD, (^), (pD, dn), (|i6i, 

and the additional 6 terms containing an e^^^^ given by Eqs. (^Bj), ( |111| ), ( |112| ), ( |125p , 
, ( P35| ). For completeness, these terms are given explicitly in Table |VII1| in Appendix 



Eqs. 



1), and (m, had 



0. Recall that three additional equation of motion terms, 
previously been eliminated using the trace relations. 

It is important to note that this argument must be applied with extreme care if one works 
in the other direction. Thus if one starts with a particular Cq generated from some model 

C2) 

which contains terms proportional to O^qj^,^ and tries to put it into the general form we 
have derived it will be necessary to keep track of the changes in the coefficients generated by 
the successive transformations, particularly the second order correction to £2 which would 
be missed by just dropping the C^om terms. This is explained in more detail in 
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B. Reorganization 



In the course of the derivations of the preceding sections we devoted a great deal of effort 
to getting all of the terms of the Lagrangian in a systematic way, without any thought at all 
as to which terms would be of most practical importance. It is clear however that there will 
be some terms which will be of immediate importance for simple processes. In fact there 
have been already recent calculations, e.g. Ref. ||3^, which in the absence of the general Cq, 
have included a few ad hoc terms, motivated by the need to cancel infinities arising from 
the loops involving £4. On the other hand some terms in Cq contribute only to processes 
which are so complicated that they probably will not be of practical interest for a long time. 
For example, naively the factor [D^U]- goes like d^(j) in leading order and so it would seem 
that a term like that of Eq. ( |5^ ) would contribute at tree level only to a process involving 
six boson fields, which is probably not of much immediate interest. 

One concludes from the preceding discussion that it would be useful to organize the terms 
in the Lagrangian in such a way as to separate out those contributing to simple processes. 
Such an organization is more subtle than it might seem however. Consider for example the 
simplest factor appearing in many terms of the Lagrangian, [D^U]^. If we start with Eq. 
(p!5|) for U and expand in powers of (j) we get symbolically [/ ~ 1 + i(f> + 0{(t)'^), where we 
have absorbed the Fq into (j) for the purposes of this section. Now using the first line of Eq. 
(p!8|) we expand [D^iU]_ as 

[D,U]^ ~ id,4> - i{R, - L^) + [L^, 4>] + (263) 

If = = 0, corresponding to pure QCD with no external fields, then [D^U]^ ~ d^(j) 
and a term like Eq. (|5^) which involves six D^'s does contribute only to a process with six 
bosons in accordance with our naive expectation. In general however when 7^ this 
term could contribute to anything with n bosons and 6 — n external fields, where < ?7, < 6. 
Thus for the general case it doesn't appear possible to sort the terms in a useful way. 

However if we limit the external fields to the electromagnetic field, which is a useful phys- 
ical situation, then we can take R^ = = —eA^Q, where is the electromagnetic field, 
e > is the electric charge, and Q is the diagonal quark charge matrix 3Q = diag(2, —1, — 1). 
This allows us to classify the various terms according to the number of boson and external 
electromagnetic fields required in the process in order that the term produce a non zero 
contribution in this limit. The explicit term of course remains perfectly general, but the 
sorting into various groups depends on this special assumption. 

With the assumption of only external electromagnetic fields, i.e. i?^ = = —eA^Q, 
and with x = the various building blocks contribute to order and lower as: 





~ 0, A(f), 




~ 0, A(f), A^(j), 


[xU 




[D.xU 


~ x,x(P,Ax,Ax<P, 


M- 


~ x4>, 


[D,X]- 


~ x4>, ^x4>, 
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_ ~ A,A(f), 








_ ~ A.Acj), 




~ A(j), 











(264) 

Note that for the electromagnetic case and in addition the usual choice of x = '2BqM, 
with M the diagonal quark mass matrix, M = diag(m„, m^, m^), the quantities [-D^x]± and 
[D^D^x]± actually vanish since [Q,x] = implies D^x ~^ d^x- We keep them, however, to 
preserve a bit more generality. 

There is another useful simplification we can use also. Under parity, using the results of 
Sec. |II E 5| , we have in effect for each term U U'^ ^ x ^ x\ and F^^, 

with an extra minus sign for the terms. This means that for the electromagnetic or 

pure QCD case, with = and Fj^'^ = F^" , and with x = x\ we can use the fact that 
[/ — s> f/"!" is equivalent io (j) ^ —cj) to show that parity invariance implies that terms without 
an Cafi-yS will have only even powers of while those with an Cq/j^^ will have odd powers 

US. 



We have used both of these simplifications to group the terms in the final result of Tables 
P|- |VII| according to the smallest number of ^'s and A^'s which a term can have. There can 



be more than the minimum number of course, if one goes past leading order in the expansion 
in powers of (f). Also it is possible that there may be accidental cancellations which make 
the leading term vanish so that a given structure may in fact have more (/)'s or A's than 
indicated in the tables. In a few cases the leading behavior of the structure contains no 
0's. These terms have been included with the terms having two 0's, since such non leading 
contributions would seem most relevant for practical calculations. 

To see how this sorting works, suppose we are interested in the contact terms at order 
contributing to the process 7 + 7 — > vr + vr. We thus need two 0's and two A^'s in tree level. 



Thus none of the terms in Tables |I| - [VI I| contribute as they require an odd number or too 
many (;/)'s. Of the terms in Table |I| the first group has no A^'s and the last group has too 
many, so neither will contribute. Thus we need consider only the middle three groups, and 
even some of those terms will vanish because as noted above [£)^x]± and [D^Dyx\± vanish. 

If one wishes to consider a general external interaction which has 7^ then there 
seems to be no substitute for considering each term in detail to determine which will con- 
tribute. We emphasize again that the results listed in Tables || - [VI I| are completely general 



and appropriate also for general external interactions. Only the classification into groups 
depends on the assumption of just the electromagnetic interaction. 

C. Simplifications 

The notation used so far was developed to simplify the derivations. It has the advantage 
of leading to a Lagrangian which depends on a relatively small number of building blocks 
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which have well defined and simple transformation properties under parity and charge con- 
jugation and hermitian conjugation. 

For the purposes of calculation however there are a few simplifications which will be 
collected here for reference, though most have been mentioned earlier. Some may simplify 
the evaluations in specific cases. 

The building block [D^U]- can be written as [D^U]- = {D^U)W = -U{D^U)^. Under 
the trace the t/'s and t/^'s for the most part commute through and collapse to unity so that 
the effect is to convert a string of [D^[/]_'s to a string of (D^f7)'s and (D^f7)''''s. 

The factors and [H^'^]^ can be expressed in terms of the original and F^'^ via 

[G^^]^ = G^'^'U^ = UG^""^ = {F^^ + UF^^U^), (265) 

[H^"']+ = H^'^U^ = UH^""^ = {F^" - UF^^U^). (266) 

For purely electromagnetic external gauge fields, with proportional to the (diagonal) 
quark charge matrix, or for pure QCD with no external gauge fields, and for the usual choice 
for X as a diagonal quark mass matrix the covariant derivative Dfj,x — >■ — > 0. Thus in 
this situation the terms containing [-D^x]± or [DnD,yx]± ^ill vanish. 



D. Final Results 



We have collected our final results for Cq, the complete Lagrangian to order p^, in Tables 
m - [VI I|, ordering the various terms according to the scheme described in Sec. [IV B| above. 



Each term has been written in such a way that it is chirally invariant and for real coefficients 
is hermitian and invariant under parity and charge conjugation. There are a total of 32 terms 
of odd intrinsic parity, corresponding to the coefficients Ai and involving an ea/3'ys- There 
are 111 terms of even intrinsic parity, corresponding to the coefficients Bi. In the course 
of the derivation we obtained 23 independent equation of motion terms proportional to 
the operator C^oa/ of Eq. ( [260[ ). For completeness these have been listed in Table [VIII 



in Appendix though we assume that for the simplest Lagrangian these will have been 
transformed away via an appropriate field transformation . 



Also in the course of the derivation we used trace relations to express 18 structures 
originally obtained in terms of others. Table |^ of Appendix ^ shows which of the original 



terms are related and which we chose to eliminate. Likewise epsilon relations were used to 
eliminate 16 dependent structures from the original set of odd intrinsic parity terms. Table 
^ of Appendix indicates which of the original equations were eliminated. Both of these 
tables, though not strictly necessary for the final results, should make it easier to compare 
our work with that of others. 

We have tried to start with all possible structures and to eliminate those which are not 
independent and to extract from the remaining terms as many equation of motion terms 
as possible. The procedure for doing this requires several tricks, i.e. the trace and epsilon 
relations, and depends in some cases on the way in which one writes the various terms. We 
cannot prove in a general and rigorous way that the resulting terms are all independent. 
Thus the reader should be aware that it is possible that there may be additional tricks which 
have been missed which could be used to express some of the structures in terms of others 
and thus reduce the number of independent terms. 
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E. Comparison with other results 



The set of independent structures of 0{p ) in the odd intrinsic parity sector has aheady 



been discussed by Issler |^ and by Akhoury and Alfakih |]36|. However, the number of 
independent terms we find does not agree with either of the above references, which mutually 
disagree with each other. In the following, we will try to locate potential sources of this 
discrepancy. A direct term-by-term comparison is made difficult by the fact that, in general, 
different conventions and, furthermore, different basic building blocks are used. 

Let us start with Ref. which quotes 49 independent terms. This number is close to 
our starting number of 54 terms. However, we made use of the epsilon relations to eliminate 
16 terms. It appears that these relations were not used in any derivation prior to the work 
by Akhoury and Alfakih (see for example Thus one has to conclude that, in 



general, too many terms were found which, in fact, are not independent. Furthermore, 
there is no reference in Bl] to the use of the equation of motion or field transformations 
to eliminate terms. For example, of the 4 terms proportional to kf"* — kf^ in Ref. [0], the 
epsilon relation can be used to eliminate say the structure proportional to kf'\ and from 
the remaining 3 terms the ones proportional to k^^^ and ^2°'' can be related using a field 
transformation, resulting in only two, instead of four, independent structures. We made use 
of 6 field transformations which reduces our final number to 32 as compared with 49 in Ref. 

However, even after taking these two observations into consideration, there remain some 
discrepancies. It appears that in Ref. not all independent orderings of operators under 
the trace have been taken into account. As an example, in our opinion, there should be 
another term similar to the structure proportional to ^[2'' involving a different contraction 
of indices. 

Finally, it appears that the set of terms includes structures which can be related using the 
total derivative argument resulting in a reduction of the number of terms. To be specific, let 
us consider as an example the structure proportional to ky^. It is straightforward but tedious 

to show that up to a total derivative it is related to the terms proportional to k^^ , k^^ , k^J 

(2) 

and kii . For that purpose one has to take the covariant derivative off of the field strength 
tensor in the term proportional to k^^^ use the total derivative argument as outlined in 

Sec. ^nn\ 

A comparison with the work of Akhoury and Alfakih |36| turns out to be more difficult 
as their choice of the building blocks is very different from ours. The final number quoted 
in Ref. [^] is 30 where 5 terms have been eliminated using the equation of motion. This 



has to be compared with our 32 terms using 6 equation of motion terms. 

Even though the use of epsilon relations was first proposed in Ref. it seems that 
their set still contains structures which are not independent as a consequence of such re- 
lations. To give an example, the terms proportional to Wu — W13 can be interpreted to 
originate from a tensor Qa/s^iupa which is antisymmetric in the index pairs (a, (/i, u), and 
(p, cr), respectively. Without the epsilon relation, one would naively expect 3 independent 
contractions from such a tensor, of which only one independent term remains after use of 
the epsilon relation. In a similar fashion one can show that of the four terms proportional 
to Wj — WiQ only three are independent. Finally, the very first structure proportional to Wi 
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vanishes identically (see our Eq. ( P58| )). 

On the other hand, it appears that there are terms missing in Ref. [^]. To be specific, 
there exists an additional independent contraction of indices for the structure of the type 
proportional to Wq. Furthermore, Ref. does not contain any terms involving covariant 
derivatives of x of which we find three independent terms. Finally, note that in Ref. the 
equation of motion is used so as to eliminate terms which are proportional to instead 
of the structures proportional to [Df^D^U]-. On the other hand this means that in Ref. 
terms proportional to should have been kept. In fact, we find six such terms (see 

Table |V11J| ) whereas Ref. quotes only three. 

Thus it appears that in both of the previous cases where a systematic study of the 
odd intrinsic parity terms was made there are terms in the resulting sets which are not 
independent and terms which have been missed. 



V. SUMMARY 

In the preceding sections we have developed the complete chirally invariant Lagrangian 
Cq for the meson sector to order p^. This is intended to be an extension of the order p'^ 
Lagrangian £4 of Gasser and Leutwyler which has become the standard in chiral perturbation 
theory and has been used in many applications. Such an extension is important at this time 
because we are beginning to see two loop calculations of processes for which the leading 
contributions are order p'^. Such calculations generate some p^ contributions, but the full 
Cq is needed to produce a consistent result. 

Throughout we have emphasized a careful and pedagogical development of the steps 
leading to the full Lagrangian, since we feel that it is only via such an approach that the 
reader can be confident that the extremely complicated final result is complete and correct. 
To do this we have first outlined a hierarchical strategy which allows us to eliminate terms 
in favor of ones lower in the hierarchy. We then discussed a number of general results which 
allowed us to simplify and reduce the number of terms. After imposing parity and charge 
conjugation invariance we could obtain a set of general structures at each level, which could 
then be evaluated to find the set of possible terms. Trace relations and epsilon relations 
were then used to eliminate terms which were not independent and we described how field 
transformations could eliminate those terms proportional to the lowest order equation of 
motion. 

The resulting set of terms was then sorted, for the usual QCD plus electromagnetic case, 
according to the minimum number of boson and electromagnetic fields appearing. The final 
result for Cq is given in Tables ^ - |V11| . It consists of 111 terms in the even intrinsic parity 
sector and 32 terms in the odd intrinsic parity sector. 

To our knowledge there have been no prior systematic studies of the even intrinsic parity 
sector to this order, though isolated terms have been used in a variety of calculations. In 
the odd sector however there have been two previous analyses [04] , ^ , which disagree in the 



number of terms with our result and with each other. We have shown that in each of these 
previous cases, terms have been missed and terms which are not independent have been 
included. 
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It is clear that the coefficients of all of these terms will never be evaluated from experi- 
ment. However a much smaller subset actually contributes to most simple processes, and it 
may be possible to get information about some of them. 

In any case we hope that our derivation of the complete and most general Lagrangian 
will stimulate systematic chiral perturbation theory studies of processes at this order. 
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APPENDIX A: RELATIONS BETWEEN TRACES 



Following the lecture notes of Coleman |ED[, we will derive relations between traces of 
3x3 matrices. Let A be any complex 3x3 matrix with eigenvalues ai, 02, and 03 (possibly 
complex and identical). The solution of the characteristic equation is then equivalent to 



{A - aiI){A - a2l){A - Usl) = 0, 



(Al) 



where / is the 3x3 identity matrix. As A is similar to a matrix B of the form 

5, 5' = Oorl, 





1 


ax 


b 





A = TBT-\ B = 









b' 




\ 








as 



(A2) 



one finds 



Tr{A) 
Tr{J^) -- 
det{A) 



= Tr{B) 
Tr{B^) 
-- det{B) 



ai + 02 + as, 
al + al + al, 
a\a2a^ . 



These relations may be used to rewrite Eq. (|A1|) as 

A^ - Tr{A)A'' + - {iTr{A)f - Tr{A^)) A - det{A) = 0. 



(A3) 



(A4) 



A ffist important observation is made by taking the trace of Eq. (|A^), namely that the 
determinant of a matrix can be expressed in terms of traces. This is the justification for 
not considering determinants as separate building blocks in the construction of the chiral 
Lagrangian Taking the trace of Eq. (|A4|) one eliminates the determinant to obtain (see 
also Eq. (80) of Ref. R]) 



A^ - Tr{A)A^ + ^{Tr{A)yA - ^Tr{A^)A 
-^Tr{A') + ^Tr{A')Tr{A) - ^{Tr{A) f 



(A5) 
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Starting from Eq. ( |A5|) we will derive various trace relations for traces involving between 
four and six 3x3 matrices. 

Multiplying Eq. ( [A5| ) by A and taking the trace results in 

Tr{A^) - -Tr{A'^)Tr{A) - -{Tr{A^))'' + Tr{A^){Tr{A)y - -{Tr{A))^ = 0. (A6) 
3 2 6 



Inserting A = XiAi + X2A2 + A3A3 + X^A^ into Eq. ( |A6| ) and comparing the coefficients of 
A1A2A3A4 one finds (see also Eq. (81) of Ref. [^ ) 



J2 Tr{A,A2AsA^) - Tr{A,A2A,)Tr{A,) - Tr{A,A2)Tr{A,A^) 

Sperm. Sperm. 3 perm. 

+ Y Tr{AiA2)TriA3)Tr{A^) - Tr{A^)Tr{A2)Tr{As)Tr{Ai) = 0. 

Sperm. 

(A7) 

In the following, we list special cases of Eq. ( |A7| ) which we used to relate different terms 
and thus eliminate redundant structures in the chiral Lagrangian: 

J2 Tr{A,A2AsA^) - Tr{AiA2As + A,A,A2)Tr{A,) 

Sperm. 

- Y Tr{A^A2)Tr{AsA,) = 0, (A8) 

3 perm. 

for Tr{Ai) = Tr{A2) = TrlA^) = 0, and A4 arbitrary; 

Y Tr{AiA2A3A^) - Yl Tr{A,A2)Tr{AsA^) = 0, (A9) 

Sperm. 3perm. 

for Tr{Ai) = 0; 

2Tr{{A^B + ABA + BA^)C) - 2Tr{A^B)Tr{C) 

-Tr{A^)Tr{BC) - 2Tr{AB)Tr{AC) = 0, (AlO) 

for Tr{A) = Tr{B) = 0, and C arbitrary; 

2Tr{{A^B + ABA + BA'^)C) - Tr{A^)Tr{BC) - 2Tr{AB)Tr{AC) = 0, (All) 

for Tr{A) = Tr{B) = Tr{C) = 0; 

ATriA^B^) + 2Tr{ABAB) - ATr{A^B)Tr{B) 
-Tr{A^)Tr{B^) - 2{Tr{AB)f + Tr{A^){Tr{B)f = 0, (A12) 

for Tr{A) = 0, and B arbitrary, and 

ATriA^B^) + 2Tr{ABAB) - Tr{A'^)Tr{B^) - 2{Tr{AB)f = 0, (A13) 

for Tr{A) = Tr{B) = 0. 
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The last relation, Eq. ( |A13|) , was already used by Gasser and Leutwyler in the construc- 
tion of the p'^ Lagrangian Note, however, that the matrices A and B do not have to be 
hermitian for Eq. ( |A13| ) to hold, as is sometimes stated in the literature. Furthermore, Eq. 
( |A9| ) is the result obtained in Ref. . 

Next we multiply Eq. ( [X5| ) by A"^, take the trace, and rewrite Tr{A'^) using Eq. (^) to 
obtain 

Tr(^5) - ^Tr{A^)Tr{A'') - ^Tr{A'^){Tr{A)f + ^Tr{A^){Tr{A)f - ^{Tr{A)f = 0. 

(A14) 

Inserting A = X]^=i \Ai into Eq. (|A14|) one finds 

^ TriA^A^AsA^A,) - ^ Tr{A,A2A,)Tr{A,A,) 

24 perm. 20 perm. 

- Yl Tr{A^A2A,)Tr{A,)Tr{A,) + 2 Yl Tr{A^A2)Tr{A;)Tr{A,)Tr{A,) 

20perm. lOperm. 

-ATr{Ai)Tr{A2)Tr{A3)Tr{Ai)Tr{A5) = 0. 

(A15) 

We applied the following special case of Eq. ( [A15| ) 

Y Tr{AABBC) - -Tr{A^C)Tr{B^) - -Tr{B^C)Tr{A^) 

6 perm ^ ^ 

-Tr{ABC)Tr{AB) - Tr{ACB)Tr{AB) 
-Tr{A^B)Tr{BC) - Tr{AB'^)Tr{AC) = 0, (A16) 

for Tr{A) = Tr{B) = 0, and C arbitrary. 

For our final application we can restrict ourselves to Tr{A) = 0. Multiplying Eq. (|A5|) 
by A'^, and taking the trace one obtains 

Tr(A^) - -Tr{A^)Tr{A^) - -{Tr{A^)y = 0, (A17) 
2 3 

for Tr{A) = 0. Inserting A = Y.Li ^i^i with Tr{Ai) = into Eq. ( [AT7| ) yields 
Y Tr{A,A2A,A,A,Ae)-'^ ^ TriA,A,A3A,)Tr{A,A,) 

120 perm. 90 perm. 



- J2 TriA,A2A,)Tr{A,A,A,) = 0, (A18) 

40 perm. 

for Tr{Ai) = 0. As a special application we insert Ai = Pa, A2 = P", A3 = Pp, A4 = P^, 
^5 = P^, and Aq = P^ into Eq. (|Al8|) to obtain 

2rr((P • Pf) + 3Tr(P ■ PP^P ■ PP^) + 6rr(P ■ PP^P^P'^P") 
+3Tr{P^P^P^'PpP''PP) + TriP^P^PpP^^P^Pf) 

-Tr{P ■ P)Tr{{P ■ Pf) - ^Tr{P ■ P)Tr{PpP^P^'P'') 

-ATr{P^P^)Tr{P ■ PP'^P") - 2Tr{PpP^)Tr{PPPpP'' PP) 
-3Tr{P ■ PP^)Tr{P ■ PP^) - Tr{P^,P^Pp)Tr{PPP'' Pf) 

-Tr{P^P^Pp)Tr{P''PPP'') = 0, (A19) 
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for rr(P„) = 0. 

In Table |X| we summarize how we applied the trace relations. It contains the equation 
numbers of structures which are not independent due to trace relations as well as the specific 
trace relation which connects them. Furthermore, we list which structure we have chosen to 
eliminate. 



APPENDIX B: EQUATION OF MOTION TERMS 

For the purposes of completeness we list in Table [VTTT| in the same form as our final results 
all of the terms in the original Lagrangian which are proportional to the factor [Df^V^U]-, 
which has been replaced by O^^qj^j. As detailed in Sec. |IV A| these terms can simply be 
dropped from the most general form. 



APPENDIX C: EPSILON RELATIONS 

We list in Table ^ the equation numbers of the original structures we derived which are 



related by the epsilon relations of Sec. II E 7 and the ones which we chose to eliminate. 
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TABLES 



TABLE L Transformation properties under the group (G), charge conjugation (C) and parity 
(P). The expressions for adjoint matrices are trivially obtained by taking the hermitian conjugate 
of each entry. In the parity transformed expression it is understood that the argument is (— x, t) 
and that partial derivatives 5^ act with respect to x and not with respect to the argument of the 
corresponding function. 



operator 


G 




C 




P 


U 


VrUVI 












VrDx-,...Dx^UVI 


{Dx, 


■■■Dx.Uf 




. ..L>^"C/)t 


X 


VrxVI 




X^ 




X^ 


Dx,...Dx„x 


VrDx,...Dx„xVI 


(Dx, 


■■■Dx^xf 




...D^"X)^ 




VrR^vI + iVRd^v^ 




~% 
















Ri^ 


pR 










Fl" 


pL 










F^R 




VrG^uVI 




r 






^Ai • • • Dx„G^u 


VrDx,...Dx„G^,vI 


-{Dx, 


■ ■ ■ F)xr,G^yY 


(D^i . 


B^^G^^)'^ 




VrH^^^vI 










. . . Dx„H^y 


VrDx,...Dx„H^,vI 




■ ■ Dx^H^y) 




P)Xn, J^IJ.vy{ 



TABLE IL Terms in the Lagrangian with 2 or more (/>'s, sorted according to the minimum 
number of electromagnetic fields A^. The formulas are general, but the sorting order and number of 
0's depends on the specific assumption of external electromagnetic fields only, i.e. on the assumption 



Rfj. = = —eAf^Q, and on x = as detailed in Sec. IV B . The quantities [D^x]± and [Dfj_D^x]± 
vanish if one in addition makes the usual choice x = 2BqM, but have been kept here for generality. 
Terms with leading behavior independent of 4> have been included with the two (j) terms. The double 
covariant derivative D^D^A is assumed to be symmetric in its indices in accord with Eq. (p^). 



Terms in the Lagrangian with 2 or more 0's and no A,, 's possible 



+BlTr([x] + [x] + [x] + ) ^ 



+B2Tr{[x]-[x]-\x] + ) m 



+BiTT([x] + )Tr([x] + [x]+) Mt 



-BiTr{[x]+)Tr([x]-[x]-) Ml 



-BsTr{[x]-)Tr{[x]-[x] + ) m 



+BeTT{[x] + )Tr([x] + )Tr{[x] + ) m 



+BTTr([x] + )Tr([x]-)Tr([x]-) 244 
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Terms in the Lagrangian with 2 or more 0's, with an possible but not necessary 



'U]-)) 



+BaTr ( 
+Bc,Tr ( 
+BioTr 
+BiiTr 
+Bi2Tr 
+Bi3Tr 
+Bi4Tr 
+Bi5Tr 
+Bi6Tr 
+BnTr 
+BisTr 
+BigTr 
+B2oTr 
+B2iTr 



Tr{[D^x]+) 



DcDfiU]-{[D°'U]- 
DcDpU]-lD"U]- 
ID^X] + [D>^X] + ) 

[D^x] + )Tr{[Di-x] + ) 
[D^x]-)Tr{[D'-x]-) 
ID^U]-(ID'^X] + [X]- - Ix] 
lD^U]-lD'^x] + )Tr{lx]-) 
[D^U]-lD'^x]-)Trilx] + ) 
[D^U]-lD'^U]-lx] + [x] + ) 
lD^U]-lD'^U]-[x] + )Tr{lx]+) 
[D^r/]_[DMC/]_)Tr([x] + [x]+) 
[D^U]-lx] + )Tr{lD'^U].lx] + ) 
[D^UU[D^U]-)Tr([x] + )Tr{[x] 



[Df'x]- + [x]-[D''x]+-[D''x]- 
■TrilD^U]-lx]+)Tr{lD'^x]-) 
■Tr{[D^U].lx]-)Tr{lD'^x] + ) 



M+)) 



L03 

m 
m 
m 
m 

L71 



Terms in the Lagrangian with 2 or more 0's and at least 1 Afj 



+iB22Tr ([D^D'3C/]_([DaC/]-[D,G'"']+ - lD,Gf'-^] + [D^U]-)] 
+iB23Tr llD^D'^U]-([D,U]-lDpGn+ " [DuGn + iD.U]-)) 
+iB24Tr{[D^U]-[D.U]-{[x] + [Gn+ + [Gn + [x] + )) 
+iB25Tr ( [D^ U] - Ix] + [D, U] _ lGn + ) 
+iB26Tr {[D^U]-lD,U]-[Gn + ) Tr {[x] + ) 

+iB27Tr{[D^U]-{lD,x]-[Gn+ " [x]-[D.Gn+ " [Gn + [DuX]- + [D^Gn + lx]-)) 
+iB28Tr{[D^U]-{lD,x] + [Hn+ " [x] + [D.Hn+ + [Hn+[D.x]+ ~ [D,Hn + [x] + )) 
+iB29Tr {lD^U]^lHn + )Tr {[D.xh) - Tr i[D^U]-[D,Hn + ) Tr + ) 



Eg. 
~123 
L24 

m 
m. 

m 
m 

175 



Terms in the Lagrangian with 2 or more 0's and at least 2 A,, 



Eq. 

m 

L9G 
L91 
192 
193 
195 
197 
198 
199 
20G 
Ml 
102 
122 
123 
124 
»2 

M 

m 

25C 



+B3oTr 
+B3iTr 
+B32Tr 
+B33Tr 
+B3iTr 
+B35Tr 
+B36Tr 



\\G^+ 



lD„D-yGp^] + 

[D^D-'H^^] + 
.[Di'U]-[G°'P] 

.[Df^m-lGii^] 



-[Gafi]- 



+B37Tr( [D^ U]-lDi^U]-)Tr( [G"''] + [G^n] + ) 
+B38Tr ([D^U]-[G°'I^] + ) Tr ([Dt'UUlG^^U 
+B39Tr l[D^U]^[Dl^U]-) Tr ([G"^] + [G^^] + 
+B4oTr ( [Dc U] _ [G-''] + ) Tr ( [D" U] _ [G^^] + ) 
+B4iTr ([D°'U]- [G'3^] +) Tr ( [Dp U] _ [G,,^] + ) 
+Bi2Tr ([Df,C/]_([D^G"'3] + [//„^]+ - [G^^UlDi^H^I^U 
+Bi3Tr {\D^U]-([D^G°^l']+[H^p]+ - [G"P\ + [D^ H-,p] + 
+B44Tr l[D''U]-{[D„G°'l^]+[H^p]+ - [G^l^j + iD^H^p]^ 
+B45Tr ([Z)„D^C/]_{[G-'3]+[H^^]+ - [H^^] + [G"'3] + )) 
+B46Tr([x] + [Hn + [HM 
+B47Tr{lx] + lGn + [G^.] + ) 
+B48Tr{lx]-ilHn + [G^.]+ - [G^,] + [Hn + )) 
+Bi<>Tr ilxU) Tr(lHn + [H^..] + ) 
+BsoTr{lx] + )Tr{[Gn + [GM 



^ [H^ls] + [D^^G°'^U 
[H^I3] + [D'<G°'''] + 
[H^p] + [D^G'^P] + 



[DMH-'3] + [G„^] + )) 



.[G°'i^] 
.[G"l^] 



Terms in the Lagrangian with 2 or more 0's and 3 A^'s 



riBsiTr {[Gt""] 
-jB52Tr{[Gf"'] 



]+) 
']+) 
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TABLE III. Terms in the Lagrangian with 4 or more (j)^s, sorted according to the minimum 
number of electromagnetic fields Af^. The formulas are general, but the sorting order and number 
of (/)'s depends on the specific assumptions described in Sec. IV B| and the caption of Table ||. The 



double covariant derivative DnDyA is assumed to be symmetric in its indices in accord with Eq. 



Terms in the Lagrangian with 4 or more 0's, 



+B53Tr ( [D^D, U] .[Di'D-^U]- [D^ U] _ [D" U] _ ) 
+B54Tr ilD^D, U] _ [D'^D'^ [/] _ [D'^ (7] _ [D„C/] _ ) 
+B55 Tr {[D^D, U] .[Df^D^U]- [Do, U]-[D-'U]-) 
+Bs6 Tr {[D^ t/] _ [D'^D" i7] - ) Tr ( [Do, U]-ID"U]-) 
+B57Tr {[D^D, U]-[Di'D°'U]-)Tr{[D-'U]- [D<, U] _ ) 
+B58Tr ilD^D, U] _ [Do, [/] _ ) Tr ( [D^D'' !7] _ f/] _ ) 
+B59Tr {[Df,D^U]-[D''U]-)Tr {[Dt' D"U]-[DcU]-) 
+B6oTr ( lDf,D^ U] _ [Dc [/] _ ) Tr ( [D^D" t/] _ [D" C/] _ ) 
+BeiTr ( [Dc^Dp U] _ [D" U] _ [x] - [D^U]-) 
+B62Tr} [x] + [Do, U]-[D°'U]- [D^ C/] _ [D^ (/] _ ' 
+B63 Tr ( [x] + [Da U] - C/] _ [D/^ C/] _ [D" (/] _ 
+B64Tr ( [D,Z)^i7]_[D«;7]_) Tr ([D/'C/l-M 
+B65Tr hDo,D^U]-[x]-) Tr {[D"U]-ID^U] 
+Be6Tr I [Do,D0 U] _ [D« [7] _ [D^C/] _ ) Tr ([x] - ) 
+B67Tr ( [x] + ) Tr ( [D„ U]-[D"U]- [D^ C/] _ [/] _ ^ 
+B68Tr ( [x] + [D^ U]-)Tr{[D"U]- [D^ U]-[Dl^U]-\ 
+B69Tr ( [x] + [D^ [/] - [D" (7] - ) Tr ( [D^ C/] _ [D/S [/] _ ) 
+B7oTr {lx\ + [Dc,U]-lDfsU]-) Tr ([£)-[/]_ [£)"[/]_) 
+B71 Tr ( [x] + ) Tr ( [D^ U]-ID°'U]-) Tr { [D^ U] _ U] _ ) 
+B72 Tr ( [x] + ) Tr ( [D„ U] -[DpU]-) Tr ([D''U].[Dl^U]J 

+B73Tri[D^U]-[D'^U]-[x]-[x]-) 
+Br4Tr ( [D^ U] _ [D'^ U] - [x] - ) Tr ( [x] - ) 
+B75Tr ( [D^ U]-lDi^U].)Tr{ [x] - [x] - ) 
+B76Tr ( [D^ U] _ [x] - ) Tr ( [D" C/] _ [x] - ) g 
+B77Tr ( [D^ (7] _ [DM {/] _ ) Tr ( [x] - ) Tr { [x] - ) P72 



Terms in the Lagrangian with 4 or more (p's and at least 1 Eq. 
+iB78Tr ([D^D;3C/]-([D''C/]-[D,l/]_[H'''']+ + [/7M''] + [D,{7]_[D'3!7]-)) 

+iBrgTr ([D^D0;7]_([D,r/]_[D'3{/]_[H'^-]+ + [/fM''] + [D/5(7]_[D,!7]_)) L2G 

+iB8oTr([GM''] + ([D^t/]-[D^i7]_[DcC7]_[D"i7]_ - [DaC/]- [D°'!7]- [D^C/]- [D^C/]-)) L2£ 

+jB8iTr([Gn + {[T'M;^]-[DaC/]-[D,i7]_[D"!7]_ - [D,C/]_ [D, [7] _ [D'^[/]_ [D^ [/]_)) g 

+iBs2Tr ( [Gf"] + [D^ C7] _ [D,, U]-ID°'U]- [D, U]-) m 

+iB83Tr ( [Gf"] + [D,, t/] _ [Df, U] _ [D^ r7] _ [D" !7] _ ) 132 

+iB84 Tr f [D^ 1/] _ [D" I/] _) Tr ( [D^ i7] -[//''''] + ) 138 

+iBs5Tr hD^,Dl3U]-[D^U]-) Tr {[DfU]-[Hf'-'] + ) 139 

+iBseTr ( [Df^DpU] _ [/f""] + ) Tr ( [D/^ (7] _ [D^ (7] _ ) m 

+iBs7Tr ( [G""] + [D, I/] _ ) Tr ( [D" C/] _ [D^ C7] _ [D, [/] _ ) IM 

+iBssTr ( [Gf"] + [Df. U] _ [D, C7] _ ) Tr { [D^ !7] _ [D" [/] _ ) L42 

+jB89Tr([Gf"'] + ([D^C7]-[Dc<i7]- - [D„i7]_ [Df,(7]_)) Tr {[D°'i7]_ [D,l/]_) 142 
+iB9oTr([D^r7]_[D,C/]_([x]-[J/n+ " [Hn+lx]-)) 

+ iBgi Tr { [D^ U] _ [Hf + ) Tr ( [D, [/] _ [x] - ) 183 
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Terms in the Lagrangian with 4 or more 0's and at least 2 A^'s Eg. 

+B^2Trl[D^,U]-[D^'U]-[H°'t'] + [H„f)] + ) ^ 

+S94Tr ([Z),C/]_[Z)^C/]_[//^^] + [/f"T] + ) m 

+B95Tr {[D^U]-[D^'U]-)Tr {[H°'P] + [H^^]+) MS 

+S96Tr r[i)^(7]_[^f°/']+) Tr {[Df'U]-[H^p] + ) m 

+BgrTr l[D^U]-[Dl^U]-) Tr ([//"^] + [//^^] + ) 21G 

+BQsTr { [D,, (7] _ [H-^] + ) Tr ( [D" [/] _ [H^^] + ) ni 

+B99Tr ( [D- C/] _ [/f^^] +) Tr ( [D^ [/] _ [/f,^] + ) 212 
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TABLE IV. Terms in the Lagrangian with 6 or more (j)^s. The formulas are general, but the 



sorting order and number of (/)'s depends on the specific assumptions described in Sec. |1VE| and 
the caption of Table ||. 

Terms in the Lagrangian with 6 or more 0's, with an possible but not necessary Eq. 



+BiooTr ( [Dc U]-[D°'U]- [Dp U] _ [D" U] _ [D^ U] _ [D'< U] _ ) 54 

+B101 Tr f [D„ U]-ID°'U]- [DpU] _ [D^ U] _ [D" U] -[D^'U]^] 55 

+Bio2Tr hDc,U]-[D°'U]-[Di3U]-[D^U]-[D'<U]^[Dl^U]^j 5C 

+B103 Tr I [Dc, U] - [Dp U] - [D^ U]-[D°'U]-[Dl^U]- [Di U]-) 57 

+Bio4Tr {[D^U]-[Dt'U]-)Tr [[Dc,U]-[D°'U]-[DpU]-[DfU]-) 3S 

+Bio5Tr ([D^U]-[D^U]-)Tr i[Di'U]-[D''U]-[Dc,U]-[D°'U]-) E 

+Bio6Tr ([D^U]-[Di'U]-[DcU]-)Tr [[DpU]-[Dl^U]-[D''U]-) TQ 

+Bio7Tr { [D^ U]-[D,U]- [D^ U]-)Tr {[D''U]-[D-'U]-[D°'U]-) 72 

+B108 Tr ( [D^ U] _ [D, U] _ [D^ U]^)Tr {[D''U]-[D°'U]-[D''U]-.) 7J 

+B109 Tr ( [D^ U] _ [D" U] -)Tr( [D^ U]-[D''U]-)Tr [lDpU]-[Df<U]-) [5 

+B110 Tr { [D^ U] _ [Di^U] _ ) Tr ( [D,, C/] _ [D^ [/] _ ) Tr ( [D- C/] _ C/] _ ) re 

+Bi 1 1 Tr { [D^ t/] _ [D, [/] _ ) Tr ( [D*^ I/] _ [D"' [/] _ ) Tr { [D" [/] _ [D^ [/] _ ) 77 



TABLE V. Terms in the Lagrangian with 1 or more (/)'s and an e^^^^, (eoi23 = +1)) sorted 
according to the minimum number of electromagnetic fields Af^. The formulas are general, but the 



sorting order and number of i;^'s depends on the specific assumptions described in Sec. IV E and 
the caption of Table 

Terms in the Lagrangian with 1 or more 0's, an ta/l-yS, and at least 1 A^i Eq. 

+AiTr [[D^^U]-{[D-x] + [G''']+ - [x] + [D"G-i']+ - [G-^^] + [D- x]+ + [T"'G^*] + [x] + )) ^^-^7^ El 



Terms in the Lagrangian with 1 or more 0's, an eg/a-yii a^nd at least 2 A^'s Eq^^ 

+iA2Tr ([D^'U]-{[D•'G°'f'] + [Gc.'<]+ - [G°'^] + [D'^gI'<] + )) e^,p^ p 
+iA3Tr ([D'^!7]_([DcG"''] + [GT*]+ - [G"f] + [DaGT"*^]+ - [D„GT"'] + [G"'3]+ + [GT"'] + [DcG°f ] + )) e^p^s ^ 
+iAATr\[x]-[G°'^] + [G-'^] + ) e^p^s »2 
+iA5Tr {[x\ + {[H'-^] + [G'<'\+ - [G^^] + [//"'3] + )) e^p^s »3 
+iA6Tr([x]-)Tr ([G"'']+ [G^'V) ^^^^^ EI 
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TABLE VI. Terms in the Lagrangian with 3 or more (p''s and an ea/s-yS, (eoi23 = +1), sorted 
according to the minimum number of electromagnetic fields Afj_. The formulas are general, but the 



sorting order and number of i;^'s depends on the specific assumptions described in Sec. IV E and 
the caption of Table ||. The double covariant derivative 0^0^ A is assumed to be symmetric in its 
indices in accord with Eq. (27). 



Terms in the Lagrangian with 3 or more 4>'s, an ea/s-yS, and at least 1 A/j 



+A-rTr ([D°'Df'U]-{[D|3U]-[D-'U]-[G^"']+ - [G'"'] + [DTf/]_[D^C/]_) 
+A8Tr l[D°'Df^U]-{[D''U]-[Df3U]-[G^"']+ - [G^'•'] + [DpU]-[D~'U]-) 
+AgTr l[D"Df3U]-[D''U]^[G^''] + [D'^U]-) e^.,s,. 
+AioTr {[D'=^DPU]-{[DpUU[D^H'''^]+ + Hi^'^] + [DpU].)\ e^^^,^ 
+AiiTr ([D"D'3C/]_([DT-C/]_[Z3^HA"']+ + P^/^^l + pTt^]-)) cc^m- 

+Ai3Tr \\D>'U\-\x\-{D''U\^{G-^^\ + ) e^.^s 
+ Ai4Tr f [D*^ C/] _ [D- [/] _ [G^^] +) Tr { [x] _ ) e^,.,s 
+A^^TrhD^^U\-[D''U]-{[x] + [H^']+ - [H-<^] + [x] + )) e^u^s 
+AieTr i[D'^U]-[H''''] + ) Tr ([D-U]-[x] + ) e^^^^s 

+AirTr l[Df^U]-{[D-x]-[H-'']+ - [x]-[D-'H-'']+ + [H-''^] + [D-x]- - [D- H-''] + [x]-.)) ^^.-yS 
+Ai8 (Tr {[Df'U]-[H^^] + ) Tr ([D-^x]-) - Tr ([D''U]-[D'^ H-<'^] + ) Tr {[x]-)) e^^^g 



L13 



L15 



L17 



126 



L27 



184 



185 



L86 



L87 



188 



L77 



L78 



Terms in the Lagrangian with 3 or more 0's, an eg/a-yii and at least 2 A^'s Eq. 

+iAigTr 1[D^U]^IDI'U]-{IG°'^] + IH^^]+ - [H^^] + [G"''] + )) e^^^s pl3 

+iA2oTr l[D''U]-[D'^U]-{[G"^] + [H^-']+ + [H^-<] + [G°'0] + )) e^,^^ 214 

+i AaiTr ( [D'^C/] _ [G""] + [D^U] „ +) 215 

+iA22Tr ([D''C/]-([D„t/]-[G°'''] + [HT*]+ - [//T*] + [G°=f ] + [£><,[/]_)) e^fi-,s iW 

+iA23Tr hDi'U]-(lD-'H°'l^] + lHc,^]+ - + + e^,p^ 225 

+iA2iTr \[D>'U\-{[DcH°'l^\ + [H'<^]+ - [H°'l^\ + [Dc.H''^]+ - [D„Ht*] + + [Ht*] + [D„H°'3] + )) 226 

+iA25Tr([x]-[H"''] + [H^*] + )ec^.,5 251 

+iA26Tr{[x]-)Tr {[H"f<] + [H^^]+) e^p^s §54 



TABLE VIL Terms in the Lagrangian with 5 or more 0's and an CafB^yS, (eoi23 = sorted 
according to the minimum number of electromagnetic fields Af^. The formulas are general, but the 



sorting order and number of (/>'s depends on the specific assumptions described in Sec. IV B and 
the caption of Table 0. The double covariant derivative D^DyA is assumed to be symmetric in its 
indices in accord with Eq. 



Terms in the Lagrangian with 5 or more 0's and an tg^-yS) with an A^ possible but not necessary Eq. 



ViA27Tr {\D^'D''U]-{[D^U]-[Dl^U]-[D-^U]-[D^U]- + [D^U]-[D-'U]- [D^* U] - [D, C/] _ )) e^ft^s 

+iA2sTr [[x\-[Dc.U\-[D0U\-[D^U\-[DsU\-) e^^T^ ^4 



Terms in the Lagrangian with 5 or more 0's, an tafl-yS, ^Liid at least 1 A^j Eq. 

+A29Trr[HM-] + {[D"C/]_[D'3(7]_[D^C/]_[D^i7]_ + [D^(7]_[D^C/]_[D'3l/]_[D"l/]_)^ e^^^p E 
+ AsoTr ( [HA"-] + { [D° C/] _ [D^ U] _ [Of' C/] _ [DW] _ + [D^ (7] _ [D/^ [/]_ [DW] _ [D" {/]_)) e^,,,^ L34 
+A3iTr ( [//'"']+ [D-C/] _ ) Tr ( [Dl^U] _ [D^ U]-[D^U]-) 6^,^/3 L44 
+A32Tr([ffn + (P"C/]-[D^C/]- + [D^C/]-[D"C/]_))Tr {[Dl^U\-[D^U]-) e^.^p \^ 
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TABLE VIII. Terms in the Lagrangian proportional to the classical equation of motion. The 
lowest order equation of motion operator O^q^ is given by Eq. ( p60|) . The double covariant 
derivative D^DiyA is assumed to be symmetric in its indices in accord with Eq. (p^), as are all 
multiple covariant derivatives. 

Terms in the Lagrangian without an edp-fS and proportional to the classical equation of motion Eq. 
+EiTr (o^^l^jlDo^D'-DpD^U]-'^ § 

+E4Tr (o'^^Uo'-^l.j) Tr {[D0U]-[D^U]-) § 
+EsTr [o'-il^^lDo^DpU]-) Tr ([D-(7]_ [D^ (7]_) ^ 



+E-rTr [ofoM[Dc.U]-) Tr {[D'-Df^U].[DpU].) g 



-EsTr[ofoj^,[DpDl'\l 

l{2) ^t(2) 



-E,Tr [O^^UPeokM\ + ^ 
-EioTr (o^^lj^jlD0U],lx]-[D^U],^ 
-EiiTr (o'-^lj^jilD^U]-lDf>x]+ + [D^x] + [DpU]-)^ 
-Ei2Tr(o'-Z,,o'^^l^,) Tr([x] + ) 



+Ei3Tr io'-^l^,j[Df,U].j Tr {[D0U].[x]-) @ 
+E^4Tr [O^eImIx]-) Tr ( [D^ [/] _ [D/^ f/] _ ) 
+Ei5Tr (o(j%j^,[D^C/]_[D''C/]_) Tr([x]-) § 
+EieTr (o^^Ij^jID^U]-^ Tr ([D''x] + ) 

+iEnTr (o^^lj^jilD^,U]^lD,U]-lHn+ + [Hn + [D.U]-[D^U]-)J 
+iEwTr (o%j^,{[D^,U]AD.Gn+ - [D.G'"1 + [Df.C7]-)) 
+iEi^Tr [0^^l^.,[D^U]-^ Tr([D,U]-[Hn + ) 
+E2oTr (ofoK'M + \x]- + M-M + )) 
-E2iTr (o%,^^[x]+\ Tr([x]-) 
-E22Tr te>5MM-) Tr{[x] + ) 



102 



118 



122 



137 



158 



159 



16C 



-E23Tr (o'-^l^{lG''l^] + [H^0]+ - [H^p] + [G"^] + )^ §31 



Terms in the Lagrangian with an ta0-^St (^0123 = +l)i and proportional to the classical equation of motion Eq. 

+iE2iTr {o%^,[D°'U]^[D^U]-[D^U]-[D^U]-^ e^p^g @ 

+E25Tr (o'^^1,^^{[DPu].[D'<U].[G^^-]+ - [Gn + [D'<U].[D0U].)^ e^^^, [n 

+E2%Tr \0%^j{Df'V]-{G''''\AD-'V]-^ ep-,^, [l2 

+E27Tr (o'-^l^ji[Dl'U]^[D-'Hn+ + [D^H^I + P'^l/]-)) e^.,^. [i 

+iE2sTr (o'^^l^jlH-^] + lH''^] + ^ e^p^s §3 

+iE29Tr (o'-^l^.j[G°'0] + [G-^'] + ) e^p^g p35 
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TABLE IX. Application of trace relations: The first column contains the relevant equation 
numbers of structures which are related by trace relations. The second column refers to the specific 
trace relation which has been applied. The third column denotes the equation number of the 
structure which we have chosen to eliminate. 



related structures 



trace relation 



structure eliminated 



,(11), i 



,(|9D,(|iD,(|3), ©,(1 



lies 



,(|l),(l9|),(|0g),(| 

(|ol),(liolD,(|io|),(|iD 

(|ig)-([l2lD,(|l3|)-(|l4g) 
Q),(S,(|167D, {|169D,(|173) 
HD, (Il70t),(|l7|) 



,(|16|),(|166 
(|9|),([194D,(|19|),(|19 
(|l9|)-O,O-(|l02D 
(|0|),(|204D,(|0|),(|0|) 
(|20|)-(|203),(H)-(|212D 



m 

m 

(P 



0) 



(1164 



m 

(|19( 



TABLE X. Application of epsilon relations: The first column contains the relevant equation 
numbers of structures which are related by epsilon relations. The second column denotes the 
equation number of the structure which we have chosen to eliminate. 



related structures 



structure eliminated 



0,(|n§,(|iT3) 

(|n|),(lii3D,(|ni),(|iii) 

(|2|),([126D,(|12|),(|12|) 
(H),(lllD,(|l3l),(|ll) 
0,(|2T3),(|22^) 



13),(12I6D, (12T8D,(|T 
0,(|226|),(|223) 
(|2|),(|229D,(|23^) 
0,(|234|) 
(|3|),(|2lD 



([13^), pel) 



(^),(|220D 



0,(|li,(|223) 



22j 

25i 
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